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What is the best way to measure baryonic acoustic 
oscillations? 
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ABSTRACT 

Oscillations in the baryon-photon fluid prior to recombination imprint different sig- 
natures on the power spectrum and correlation function of matter fluctuations. The 
measurement of these features using galaxy surveys has been proposed as means to 
determine the equation of state of the dark energy. The accuracy required to achieve 
competitive constraints demands an extremely good understanding of systematic ef- 
fects which change the baryonic acoustic oscillation (BAO) imprint. We use 50 very 
large volume N-body simulations to investigate the BAO signature in the two-point 
correlation function. The location of the BAO bump does not correspond to the sound 
horizon scale at the level of accuracy required by future measurements, even before 
any dynamical or statistical effects are considered. Careful modelling of the correlation 
function is therefore required to extract the cosmological information encoded on large 
scales. We find that the correlation function is less affected by scale dependent effects 
than the power spectrum. We show that a model for the correlation function proposed 
by Crocce & Scoccimarro (2008), based on renormahsed perturbation theory, gives an 
essentially unbiased measurement of the dark energy equation of state. This means 
that information from the large scale shape of the correlation function, in addition to 
the form of the BAO peak, can be used to provide robust constraints on cosmological 
parameters. The correlation function therefore provides a better constraint on the dis- 
tance scale ('^ 50% smaller errors with no systematic bias) than the more conservative 
approach required when using the power spectrum (i.e. which requires amplitude and 
long wavelength shape information to be discarded). 
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1 INTRODUCTION 

Prior to recombination, the ionised plasma of electrons and 
protons was coupled to the radiation in the Universe through 
the electromagnetic interaction. The pressure exerted by the 
photons worked against the gravitational collapse of pertur- 
bations in the density of baryons. This led to oscillations 
in the photon-baryon fluid. These ripples are imprinted on 
the temperature power spectrum of the cosmic microwave 
background (CMB), as detected con vincingly for the first time 
around the turn of t he millennium ( |de Bernardis et al.ll2000l : 
iHananv et al.l 120001 ). A series of acoustic peaks has now 
been measured with impressive precision by the WMAP satel- 
lite (JBennet et al.l 12003 : iHinshaw etall |2003| . |2007| , |2008| ) 
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and the Arcminute Cosmolog y Bolometer Array Receiver 
CACBAR. JReichardt et al.ll2008r ). The positions and relative 
heights of these peaks can be used to place tight con- 
straints on the values of the fundamental cosmological pa- 
rameters, particularly when the CMB measurements are com- 
bined with measurements of the galaxy power spectrum 
(lEfstathiou et al.l |2002|: iPercival et al.lj200d : ISpergel et al ' 



12003*: 'Tegmark et al."2004':'Scliak ct al. 2005l: ISanchez et~ 
1.2006; Scliak ct al. 2006; Spcrgol ct al. 2007). 

The acoustic oscillations are also imprinted on the mat- 
ter power spectrum, albeit with different phases from the 
features seen in the CMB spectrum and with a reduced am- 
plitude, due to the small fraction of the total mass in the 
Universe believed to be in the form of baryons (|Sugivama1 
1 19951 : lEisenstein fc Hulll998l : iMeiksin et aLlll999D . Recently, 
these features have attracted a great deal of interest as a 
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potential route to measuring the equation of state of the 
dark energy, ujde = foE/pDE, where Pde is the pressure of 
the dark energy and p de is its density (Blake & Glazebrook ' 
20031 : IHu fc Hainianl 2003; Lindcr 2000 ; Seo fc Eisenstein ■ 
IWang I l200^lGuzik et al.1 |2007| ; ISeo fc Eisenstein I 



2002 



20071 : ISeo et al.ll2008l '). The acoustic oscillations in the mat- 



ter power spectrum are related to the sound horizon scale. 
If we know the size of the sound horizon through measure- 
ments of the CMB temperature fluctuations, we can use this 
scale as a standard ruler. The apparent size of this ruler de- 
pends upon the parameter u;de, as this influences the angu- 
lar diameter distance out to a given redshift. 

In this paper, we focus our attention on the two-point 
correlation function which is the Fourier transform of the 
power spectrum. The series of acoustic oscillations seen 
in the power spectrum translate s into a bu mp or broad 
spike in the correlation function (JMatsubara , 2004). Such 
a feature was detected for the first time in the correlation 
function of luminous red galaxies |Eisenstein et al.l l2005l : 
lOkumura et~al] |2008| ; lEstrada et al.ll2008l ). Evidence has 
also been found for baryon oscillations in th e power spec- 
trum of galaxy clustering llCole et al.l 20051; IH iitsi 2006; 
iPadmanabhan et al.l 120071 : iPercival et al.l l20dm i. The ap- 
pearance of the acoustic oscillations in the power spectrum 
was considered in an earlier companion paper to this one by 
lAngulo et al] (|2008l ). 

Whilst the bulk of the literature on acoustic oscilla- 
tions has focused on the power spectrum, there has been 
renewed interest recently in the use of the correlation func - 
tion to extract the sound horizon scale. lAngulo et al.l (|2005l ) 
showed that the correlation function of a sample of clusters 
at 2 ~ 1, as anticipated from the Sunyaev-Ze ldovich sur- 
vey p roposed with the South Pole Telescope (jRuhl et al. I 
[2OOJ) or a red sequence photometric survey with VISTA 
(http://www.vista.ac.uk), could potentially constrain the 
sound horizon scale to an accuracy of ~ 2%, in turn fix- 
ing the dark energy equation of state to better than 10%. 

For the acoustic oscillation approach to provide compet- 
itive estimates of the value of tode, one needs to extract the 
sound horizon scale from the measured correlation function 
to sub-percent accuracy. To achieve this goal it is essential 
to quantify any systematic deviations of the scale recovered 
from the size of the sound horizon. In order to do so in a 
way that avoids the introduction of new biases or system- 
atic errors, a complete understanding of all the processes 
that shape the observed correlation fun ction is requ ir ed. 

Based on the fitting formula of ISmith et al.l (120031 ) 
and using the position of the acoustic peak in the cor- 
relation function as an estimator o f the sound horizon 
at recombination, iGuzik et al.l (|2007l ) concluded that non- 
linear evolution biases the constraints by less than 0.3%. 
lEisenstein et al.l (|2006al ') argue that any possible shift of the 
acoustic scale is unobservably small eve n at z = 0. On the 
othe r hand, the more recent analy ses of ISmith et al.l (120081 ) 
and ICrocce fc Scoccimarrd (|2008l ) have shown that both 
large volume numerical simulations and theoretical predic- 
tions based on renormalized perturbation theory indicate 
that the position of the peak in the correlation function 
shows important shifts with respect to the linear theory pre- 
diction. If unaccounted for, these shifts bias the constraints 
obtained by using BAO measurements as a standard ruler. 

In this paper we analyse some of the different tests 



that have been proposed to extract cosmological informa- 
tion from the correlation function on large scales, with par- 
ticular emphasis on the possible systematic errors that can 
be introduced in their application. We first explain the rela- 
tionship between the size of the sound horizon and the loca- 
tion and form of the peak in the correlation function in § [S] 
We next determine the accuracy of various means of com- 
puting the matter two-point correlation function according 
to linear perturbation theory for a fixed set of cosmological 
parameters (§|3]). 

Besides the linear evolution of density perturbations, 
there are other processes that will affect a measurement of 
the correlation function obtained from the new large galaxy 
redshift surveys; non-linear evolution of the density field, 
redshift-space distortions and halo bias. It is extremely im- 
portant to test our ability to model these processes in order 
to firmly assess if the precision claimed by some future ex- 
periments is attainable with our current understanding of 
these effects. In § 3] we describe the ensemble of simulations 
we used to model these effects, and the methodology im- 
plemented to estimate the two-point correlation function. 
In § [5] we describe the details of our modelling of the cor- 
relation function and its ability to reproduce the results 
from the N-body simulations with respect to non-linear evo- 
lution, redshift-space distortions and scale-dependent halo 
bias. The measurement of the imprints of the baryonic 
acoustic oscillations in the power spectrum and correlation 
functions are affected in different ways by these problems. 
In § [6] we compare the performace of these statistics to see 
which one offers the most advantages as a tool to recover 
unbiased constraints on cosmological parameters from BAO 
measurements. Finally, in §[7] we give a summary of our main 
results. 



2 WHAT IS THE RELATION BETWEEN THE 
SOUND HORIZON SCALE AND THE PEAK 
IN THE CORRELATION FUNCTION? 

The two-point correlatio n function has a bum p on scales 
in excess of 100fe~^Mpc (|Eisenstein et al.ll2005h . It is often 
assumed that the position of this bump coincides exactly 
with the sound horizon scale. In this section, we explore the 
connection between the sound horizon, the nature of the 
acoustic oscillations imprinted on the matter power spec- 
trum and the position and form of the acoustic bump in the 
correlation function. We shall demonstrate that, at the level 
of accuracy required by forthcoming measurements of dis- 
tances using BAO, the assumption that the sound horizon 
scale is equal to the position of the peak in the correlation 
function is incorrect and introduces a systematic error in 
the distance measurement which is in excess of the expected 
random errors. 

The correlation function is the Fourier transform of the 
power spectrum (see Eq. ((1} in the next section) . The acous- 
tic oscillations of the baryon-photon fluid are imprinted on 
the matter power spectrum as a series of waves, //it was the 
case that these oscillations had a fixed wavelength and am- 
plitude, then there would be a sharp feature in the correla- 
tion function at a scale centred on the wavelength of oscilla- 
tion. However, it turns out, for various physical reasons, that 
the acoustic oscillations in the power spectrum have neither 
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Figure 1. The correlation functions computed by Fourier trans- 
forming the power spectrum obtained using the EH98 formula for 
P{k) for three cases: sohd line - a fully consistent linear theory 
P{k), dot-dashed line - a P{k) with no Silk damping and with 
dominant velocity overshoot on all wavenumbers and dashed - a 
P{k) with no Silk damping but a standard velocity overshoot. 
The arrows mark the location of the peak in the acoustic bump, 
defined as the local maximum. The shaded region indicates a 2% 
error on the sound horizon scale, which is shown by the central 
solid line. 



a fixed wavelength nor a fixed amplitude. This leads to a 
more complicated relationship between the sound horizon 
scale and the position of the peak in the correlation func- 
tion than the naive assumption of equality. Furthermore, the 
size of the discrepancy between the peak location and the 
sound horizon scale depends upon the values adopted for 
the cosmological parameters, such as the matter density. 

The physical reasons behind the a ppearance of the 
acoust ic oscillations are set out clearly by lEisenstein fc HtJ 
1(1993, hereafter EH98), who discuss the form of the power 
spectrum of density fluctuations in a universe containing 
cold dark matter and baryons. There are two phenomena 
which have a direct impact on the nature of the acoustic 
oscillations. The easiest to understand is the reduction in 
the amplitude of the oscillations with increasing wavenum- 
ber. This is due to an imperfect coupling between photons 
and baryons as the end of recombination is approached; pho- 
tons diffuse out of perturbations and the remaining coupling 
with the baryons causes a Compton drag which leads to 
some baryons being pulled out of perturbations. This leads 
to a reduction in the a mplitude o f the density fluctuation 
known as Silk damping (|Silk Ill96g ). The explanation of the 
change in the wavelength of the oscillation with harmonic, 
at wavenumbers for which ks < 10, where s is the sound 
horizon, is more subtle. After recombination, the baryons 
can 'slip' past the photons. This motion of the baryons, 
called the velocity overshoot, generates new perturbations 



(jSunvaev fc Zeldovich |[l970l : iPress fc Vishniac Ill980l 'l. Dur- 
ing recombination, EH98 argue that the velocity overshoot 
is not the dominant contribution to the growth of the bary- 
onic perturbation, which leads to the "node shift" in the 
baryonic transfer function. 

The physically motivated fltting formula for the mat- 
ter power spectrum presented by EH98 can be used to show 
the impact of the phenomena described above on the form 
of the acoustic bump in the correlation function. The conse- 
quences of Silk damping and the impact of velocity overshoot 
can be readily identified in this formula and can therefore 
be switched on or off to investigate their impact on the cor- 
relation function. We consider three cases: 

(1) The fully consistent power spectrum, with Silk damping 
and a consistent treatment of velocity overshoot. 

(2) A toy model in which there is no Silk damping and 
velocity overshoot is dominant at all fc. This corresponds to 
setting fcsiik — > 00 and s = s in eqn. 21 of EH98 (s is defined 
ineqn. 22of EH98). 

(3) A second toy model with Silk damping (the appropriate 
value of fcsiik is used) but with velocity overshoot dominant 
for all fc (again, retaining s = s in eqn. 21 of EH98). 

We then Fourier transform the resulting power spectrum 
in each case to obtain the correlation function of matter 
fluctuations in linear perturbation theory. 

The results of considering the above test cases are pre- 
sented in Fig. [ij where we assume the background cosmo- 
logical parameters to be flo — 0.25, Qh ~ 0.041 and h = 0.7, 
where Ho = 100/ikms~^Mpc~^. The sound horizon scale 
for these cosmological parameters is s = 113.8/i~^Mpc (us- 
ing Eq. 6 of EH98), and is marked by the vertical line in 
Fig. [T] Using the full linear theory P(fc) (Case 1 above) 
to compute ^(r) yields a broad peak with a maximum at 
r = 111.0/i~^Mpc. The correlation function corresponding 
to Case 2 has a much sharper bump with a higher am- 
plitude than in Case 1. The maximum of this feature is 
at the true sound horizon scale. Turning on Silk damping 
(Case 3) shifts the peak to r = 113.0/i~^Mpc. The shaded 
area around the vertical line in Fig. [l] represents a 2% er- 
ror on the sound horizon measurement. This is size of the 
random error forecast fo r the ongoing Wiggl eZ experiment 
JGlazebrook et al. IIJOOtI : lAngulo et al]|2008h . The peak of 
the correlation function bump lies outside the shaded re- 
gion. Therefore approaches which advocate measuring the 
location of the bump in the correlation function will make 
significant systematic errors if this scale is mistakenly iden- 
tified as the sound horizon scale. 

The size of the discrepancy between the location of the 
acoustic bump in the correlation function and the sound 
horizon scale depends upon the cosmology. Fig. [5] shows the 
variation of the sound horizon, s (solid line), and the acous- 
tic bump location, rp (dot-dashed line) , as a function of Sim 
for a fixed value of Q^, = 0.041. It can clearly be seen that, 
although these two scales agree for f2m > 0.4, the scale rp 
is smaller than s for values of the density parameter that 
arc consistent with current observational constraints (e.g. 
iSanchcz ct al. 2006). The dotted lines in Fig. [2] indicate a 
2% error on s. For smaller values of f2m, where the baryon 
fraction is high and then the effects described above more 
important, the deviation between s and rp can be larger 
than this 2% limit. Therefore, a simple association of the 
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Figure 2. The true sound horizon, s (solid line), and the position 
of the acoustic peak in the correlation function, Tp (dot-dashed 
line), as a function of Om for a fixed value of fij, = 0.041. The 
dotted lines indicate a 2% spread in the value of the sound hori- 
zon. The dashed line shows location of the acoustic bump when 
the EH98 formalism for P{k) is replaced by a more accurate cal- 
culation made with CAMB (see text for details). 



position of the peak with the true sound horizon will lead 
to results that are biased towards high values of Q.m- This 
highlights the importance of the correct interpretation of 
the observations in order to avoid the introduction of bi- 
ases and systematic errors when using BAO observations as 
cosmological probes. 

In the next section we will see that there are dif- 
ferences between the correlation functions computed us- 
ing the EH98 fitting formula and that obtained us- 
ing CMBFAST (ISeliak fc Zaldarriaga I Il996l ) and CAMB 
(|Lewis et al.ll2000l ). The value of rp inferred using the EH98 
approximation for P{k) is biased towards smaller scales (see 
the dot-dashed line in Fig. [2| . Nevertheless, the qualitative 
description of the processes that shift the location of the 
bump in f (r) peak away from s is still correct. 

The conclusion from this simple analysis is that, even 
when using linear perturbation theory, it is wrong to assume 
that the position of the peak in the correlation function cor- 
responds exactly with the value of the sound horizon. Such 
an assumption will lead to biased constraints on cosmolog- 
ical parameters. The correct approach is to model the full 
shape of the correlation function which we pursue in the 
following sections. 



3 THE LINEAR THEORY TWO-POINT 
CORRELATION FUNCTION 

In Section 2, we demonstrated that it is not possible to use 
the position of the acoustic peak in the correlation function 
as a proxy for the sound horizon scale. In order to obtain 
useful constraints on the cosmological parameters from the 



Figure 3. The ratio of the power spectra obtained using the 
EH98 fitting formula (dot-dashed line) and CMBFAST (solid line) 
to the one obtained using CAMB. 



correlation function, it will be necessary to model the acous- 
tic bump. In this Section, we consider the form of the mat- 
ter correlation function on large scales in linear perturbation 
theory. For a realistic survey, the impact on the correlation 
function of the non-linear growth of perturbations, redshift 
space distortions and bias will have to be taken into ac- 
count. We introduce measurements of the correlation func- 
tion from N-body simulations in Section 4 and then discuss 
a full model of the correlation function in Section 5. 

The first step to a model of ^(r) is the reproduction of 
the matter correlation function in linear perturbation the- 
ory. We obtain ^(r) by Fourier transforming the linear theory 
mass power spectrum 



ar) 



f 

J — c 



A'(fc)io(fcr)dln(fc), 



(1) 



where A'^(A;) = P{k)k'^ /{2tv^) is the dimensionless power 
spectrum and jo{y) is the spherical Bessel function. 

An accurate correlation function in linear perturba- 
tion theory therefore requires an accurate calculation of the 
power spectrum. The most commonly used codes to com- 
pute the power spectr um are: i) the approximate fo rmula of 
EH98, u ) CMBFAST (ISeliak fc Zaldarriaga Ill996l ) and in) 
CAMB (JLewis et al.ll2000l ). ISeliak et al.l (|2003D tested the 
accuracy of CMBFAST against full Boltzmann codes and 
found that the relative error in the matter power spectrum 
is below 10~^. This comparison did not include CAMB al- 
though a similar accuracy is claimed for this code. The last 
release of CMBFAST was in 2003; CAMB is under continual 
development and so we adopt the view that amongst these 
three alternatives, CAMB provides the benchmark. As we 
have already mentioned in Section 2, EH98 produced a phys- 
ically motivated functional form for the power spectrum, 
which they calibrated against CMBFAST. EH98 found that 
their formalism reproduced the results of CMBFAST at the 
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Figure 4. Upper panel: The correlation functions obtained by 
Fourier transforming the power spectra computed using the EH98 
fitting formula (dot-dashed line), CAMB (solid line) and CMB- 
FAST (dashed line). Lower panel: the residuals of the correlation 
functions obtained from the EH98 and CMBFAST power spectra 
with respect to that obtained using CAMB. 



percent level, over a restricted range of wavenumbers. EH98 
provided subroutines for their fitting formulae, which are 
now in common use and are faster than running CMBFAST 
or CAMB. 

Fig. [3] shows the ratio of the power spectra obtained 
using the EH98 fitting formula (dot-dashed line) and CMB- 
FAST (solid line) to the spectrum obtained using CAMB for 
a flat cosmological model with Om — 0.237, flh = 0.041 and 
h = 0.735 (which c orresponds to the best fitting parameters 
from lSanchez et al. (2006 )). There are clear differences be- 
tween the P{k) of EH98 and CAMB on the scales relevant to 
the acoustic oscillations that can be as large as 5%. For this 
particular model, CAMB and CMBFAST show impressively 
good agreement up to fc ~ l/iMpc~^. 

The consequences for the correlation function of these 
differences in the predicted linear theory power spectrum 
can be seen in Fig. |4l The upper panel of Fig. |4] shows 
the linear theory correlation functions obtained by Fourier 
transforming the power spectra shown in Fig. [3] according 
to Eq. IT}. It can be clearly seen that while the correla- 
tion functions derived from the P{k) computed with CAMB 
and CMBFAST are in very good agreement, the use of the 
approximate P{k) of EH98 produces a correlation function 
with a quite different shape and peak position. The lower 
panel of Fig. |4] shows the residuals of the correlation func- 
tions obtained from the EH98 and CMBFAST power spectra 
from the CAMB result. For the case of the EH98 formula, 
these residuals are comparable to or larger than the variance 
expected in future galaxy surveys (see § 14. 3|) . These results 
show that the use of EH98 formula to model the shape of 
^(r) can introduce strong biases on the obt ained constraints. 
Our findings are consistent with those of ISanchez fc Cole I 



(|2008l ) who showed that the use of the EH98 fitting formula 
to model the shape of the galaxy power spectrum introduces 
changes in the recovered value of flmh of the order of one 
sigma. The EH98 formalism is excellent for providing phys- 
ical insight into the form of the power spectrum in a cold 
dark matter universe and is the simplest and quickest of the 
prescriptions listed above to use to generate large numbers of 
model spectra. However, it was never intended to supercede 
the more accurate calculations from codes like CMBFAST. 
In order to attain the level of precision demanded by forth- 
coming BAO analyses, the EH98 formalism should be re- 
placed by the more accurate calculation of CAMB when 
modelling power spectra or correlation functions. 



4 NUMERICAL MODELLING OF THE 
CORRELATION FUNCTION 

A number of effects can alter the form of the two-point 
correlation function from the linear perturbation theory 
predictions presented in the last section: the non-linear 
growth of perturba t ions, redshift space distortions and bias 
llGuzik et all l2007l: 'Smith c t al.1 l2008l: lAngulo eTall l2008l : 
ICrocce fc Scoccimarro..200& : lSmith et al.ll2008h . In this sec- 
tion we model the impact of these processes on the two-point 
correlation function using N-body simulations. The ensem- 
ble of simulations used are described in ij |4.1l The estimation 
of the correlation function in a large number of huge volume 
simulations could be prohibitively expensive without an ef- 
ficient algorithm which we describe in § 14.21 Finally, the use 
of an ensemble of simulations allows a direct estimate of the 
errors on the measured correlation function, which we set 
out in 5 1131 



4.1 The ensemble of N-body simulations 

We use an ensemble of 50 moderate resoluti on, large vol- 
ume N-body simulations called L-BASICC II (|Angulo et al.l 
|200ct). These simulations are analogous to the L-BASICC en- 
semble of simulations employed by Angulo et al. to assess the 
detectability of the acoustic oscillations in power spectrum 
measurements from future galaxy surveys. The only differ- 
ence is that the L-BASICC II runs are based on a different 
choice of the cosmological parameters. We adopt a ACDM 
cosmology consistent with current constraints from the cos- 
mic microwave background data and large scale stru cture 
measurements jSanchez et al.ll2006l : ISpergel et al.ll2007l '). We 
assume a flat ACDM cosmological model with matter den- 
sity parameter, f2m = 0.237, baryonic density parameter, 
f7b = 0.041, scalar spectral index, n^ — 0.954, a normalisa- 
tion of density fluctuations, as — 0.77, and Hubble constant, 
h = 0.735. 

Each of the L-BASICC II simulations covers a comov- 
ing cubical region of side 1340 h~^ Mpc, with the dark mat- 
ter followed using 448'' particles. This gives a particle mass 
comp arable to that e mploy ed in the Hubble Volume simu- 
lation lEvrard et al. I l|2002l ) . The equivalent Plummer soft- 
ening length in the gravitational force is e = 200/i~^kpc. 
The volume of each computational box, 2.41 /i~^ Gpc^, is 
almost twenty times that of the Millennium Simulation 
IjSpringel et al.ll2005h . and more than three times the vol- 
ume of the catalogue of luminous red galaxies from the 
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S PSS used to make the fi rst detection of the acoustic peak 
by lEisenstein et alj (|2005|). The total volume of the ensem- 
ble is 120 h~ Gpc , more than four times that of the Hub- 
ble Volume. The position and velocity of every particle are 
stored at 3 output times (z = 0.0, 0.5, 1.0). We produce a 
friends-of-friends halo catalogue at each redshift of objects 
with ten or more particles (corresponding to a mass limit of 
1.75 X 10^'' h~^ Mq). Due to their limited mass resolution, 
it is not feasible to populate these simulations with galaxies 
using semi-analytic models. 

The initial conditions for the simulations were gener- 
ated by perturbing particles from a glass-like distribution 
IJBaugh et al.lll995l ). The input power spectrum of density 
fluctuations in linea r perturbation theo ry is calculated using 
the CAMB package of iLewis et al.l (|2000l ). A different random 
seed is used for each member of the ensemb le. The simu- 
lation s were started at a redshift oi z = 63. lAngulo et al.l 
(|2003) have shown that for this choice of the starting red- 
shift, the scales relevant for the analysis of acoustic oscil- 
lations are unaffected by any transients introduced by the 
method used to generate the initial conditions. 



4.2 The practical estimation of the two-point 
correlation function 

The number of operations required to calculate the correla- 
tion function of A'^p particles by direct pair counting scales 
as Np. This is infeasible for the large number of particles in 
our simulations on the large scales considered in our analy- 
sis. The situation is further exacerbated by the fact that we 
need to repeat the calculation many times, estimating the 
correlation function for three outputs in 50 simulations. 

We follow the appro ach introdu c ed by 
iBarriga fc Gaztafiaga I (|2002l ') and lEriksen et aTl l|2004 ) to 
speed up the estimation of the correlation function. The 
first step is to construct the density field of the simulation 
on a grid of Agrid cells using the nearest grid point (NGP) 
mass assignment scheme. Using this density grid, the 
correlation function can be estimated using 



^(r) = 



A, 



^^SiSj, 



(2) 



where Si = {iii — (n))/ (n) is the density fluctuation in the 
i^^ bin of the grid and the sum extends over the Apairs cells 
separated by distances between r — Ar/2 and r — Ar/2. This 
procedure scales as Ag^j^j which is a big reduction in time 
since usually Ag^id <C Np. We use Ag^id = 240. This method 
gives an accurate estimate of the correlation function on 
scales larger than a few grid cells. In this paper, we focus on 
pair separations in excess of 60/i~^Mpc, which corresponds 
to just over 10 cells. This approach could easily be adapted 
to work on smaller scales by using different size grids to 
tabulate the density field. 

To further speed up the estimation of the correlation 
function, we first compute and store the indices of the Ancigh 
cells which contribute to a given bin of pair (cell) separa- 
tion. This list of indices can then be translated to different 
locations on the density grid to find out which cells con- 
tribute to the estimate of the correlation function in each 
bin of pair separation. In this way, no CPU time is wasted 
in re-computing cell separations. This reduces the number of 
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Figure 5. Comparison of the correlation function of one realiza- 
tion in our ensemble estimated by direct pair counting (solid line) 
with the ones obtained by counts of cubical cells with grid sizes 
of A^grid = 160 (short-long dashed line), Afg^id = 240 (dot-dashed 
line) and Afg^id = 360 (dotted line) and the ones obtained using a 
dilute sample containing a fraction of 5% (dot-long dashed line) , 
20% (long dashed line) and 40% (dashed line) of the total sample 
(panel a). The deviations from the full (direct pair count) esti- 
mate can be better appreciated in panels (b) and (c) , which show 
the difference ^grid (^) — Cfull (^) on a much expanded scale for the 
grid and for estimates derived from random samples of the dark 
matter particles respectively. 



operations from Ag^j^j to AgridAnoigh- A further speed up is 
possible as this algorithm can be naturally divided between 
processors and run in parallel. 

The assignment of the particles to the grid represents 
a smoothing of the density field and may potentially have 
an impact on the shape of the estimated correlation function 
and hence could introduce a bias into our analysis. To assess 
the effect of this smoothing, we computed the correlation 
function of one realization in our ensemble of simulations 
using a direct pair count on all scales, which is computa- 
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Figure 6. Mean correlation functions at 2 = from our ensemble of simulations in real-space (circles) and redshift-space (triangles) for 
dark matter (panel (a)) and halos in samples 1, 2 and 3 (panels (b), (c) and (d)). The error bars show the variance from t he estimates in 
the different realizations. To highlight the acoustic peak we show ^(r) X r^'^. The results in redshift space are divided by the lKaiser I (119871 ) 
boost factor (see section [5]4]l . The results for the halo correlation functions were scaled by the bias factors shown in the annotations in 
each panel (see section FS.SI I. The solid lines show the fits to the simulation data with the model for the real space dark matter correlation 
function described in section 15.11 The dotted lines in each panel show the estimates of the variance for each sample computed using 
Eq. ^. For the halo samples the variances are rescaled by the same bias factors than the measured correlation functions. 



tionally expensive. Fig. [5ja) shows the comparison of this 
estimate (solid line) with those obtained from the counts 
of cubical cells with grid sizes of A^'grid = 160 (short-long 
dashed line), A'^grid = 240 (dot-dashed line) and A^'grid = 360 
(dotted line). The deviations from the direct pair count es- 
timate, which we refer to as the "full" estimate from now 
on, are very small and can only be appreciated on a much 
expanded scale as plotted in Fig. [5jb) , which shows the dif- 
ference 5grid('") — CfuiK?") for the different values of A^grid. ft 
is clear that on using a value of A^grid = 240, the correlation 
function obtained gives an accurate description (at the level 



of a few parts in 10^) of ^fuii(r-) over the range of scales rele- 
vant to our analysis. In particular, these deviations have no 
impact on the recovered values of fc* and a (see section [S}. 

An alternative approach to make a rapid estimate of the 
correlation function is to use a diluted sample of particles 
selected at random from the simulation. Fig. [Sja) also shows 
the results obtained using a sample containing a fraction of 
5% (dot-long dashed hue), 20% (long dashed line) and 40% 
(dashed line) of the total sample. The deviations with re- 
spect to the full estimate obtained using this approach are 
larger than in the previous cases in which a grid is used. 
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Figure 7. Comparison of the correlation matrix C'or{i,j) predicted by Eq. ^ (lower triangular parts of the matrices) and the ones 
obtained from the different realizations (upper triangular parts) for the dark matter (panel (a)), and the haloes in samples 1, 2 and 3 
(panels (b), (c) and (d) respectively). In general, Eq. ^ describes accurately the elements close to the diagonal, particularly for the dark 
matter. In the more off-diagonal elements the estimation of Cor{i,j) is noisy and the agreement is not so clear. 



This can be seen more clearly in Fig. [He) which shows the 
difference ^random ('')—Cfuii(7') for the different values of A'^ran- 
Even for a value of A^'ran ~ 0.4, the deviations from the full 
measurement are larger than in the grid based measurement, 
showing that the grid based approach is the preferred one 
when dealing with large samples where the full approach is 
not applicable. More importantly the correlation function 
estimated from the diluted samples displays a systematic 
distortion from ^fuii, as it is always below the full estimate 
and shows a different shape. This systematic deviation might 
introduce a small bias when constraining cosmological pa- 
rameters. 

We apply the algorithm described above to compute the 
two-point correlation function in real and redshift space for 



the three output redshifts {z = 0, 0.5 and 1) for each of 
the L-BASIIC II simulations. Fig. EJa) shows the mean cor- 
relation function of the dark matter from the ensemble at 
z — 0. The error bars show the variance in the correlation 
function estimated from the different realizations (see sec- 
tion |4]3]). To highlight the acoustic peak we plot ^(r) x r^'^ 
ins tead of f i ^r). T he results in redshift space are divided by 
the lKaiser I ( 19871 ) "boost factor" (see section [5]4] for further 
details on the redshift space distortion). The solid lines in 
Figure |S] show the theoretical prediction for the real space 
dark matter correlation function computed as described in 
section 15.11 



We have also measured the correlation function for ha- 
los in three mass bins at each redshift. Halo sample 1 in- 
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eludes haloes with masses 1.75 x lO" < [M/H^^Mq) < 
3.5 X 10^^, halo sample 2 covers the mass range 3.5 x 10^^ < 
(M/h~^MQ) < 5.95 X 10^^, and halo sample 3 is aU objects 
with mass in excess of M > 5.95 x lO^^/i'^M©. At z = 0, 
sample 1 accounts for roughly 50% of all resolved haloes. 
The correlation functions measured for these halo samples 
are shown in Fig.[6jb), (c) and (d). The results for the differ- 
ent samples have been scaled by the respective bias factors 
h^ = 2, 78, 3.48 and 6.38 (for samples 1, 2 and 3) as described 
in the figure caption (see section 15.51 for further details on 
the distortions due to halo bias). 

4.3 The variance in ^(r) 

In this section we compare the variance in the correlation 
function estimated from the ensemble of simulations with a 
theoretical prediction. 

The starting point in the analytical estimate is the as- 
sumption that the variance in the power spectrum is that 
expected for Gaussian fluctuations with a shot-noise compo- 
nent arising from the finite num ber of objects used to trace 
the density field (Feldman. Kaiser fc Peacocklll994) : 



ap(fc) 



fH + n 



(3) 



where V is the simulation volume, and n is the mean density 
of the object s cons idered (dark matter particles or halos). 
lAngulo et al.l (12003) found good agreement between this ex- 
pression and the variance in P{k) measured from numerical 
simulations. 

The co variance of the two-point corre lation function is 
defined by (ICohn II2006I : ISmith et al]|2008l ): 



cdry) 



((c(0-e»)(e(O-^>'))) 

dkk^ 



27r2 



■Jo{kr)jo{kr')ap{k), 



(4) 
(5) 



where the last term can be replaced by Eq. ((3)1 . The variance 
in the correlation function is simply (j|(r) = C^{r,r). The 
direct application of Eq. ^ would, however, lead to a sub- 
stantial overprediction of the variance, since it ignores the 
effect of binning in pair se paration which reduces the cova ri- 
ance in the measurement (|Cohn II2006I : ISmith et al.ll2008h . 

An estimate of the correlation function in the i^^ pair 
separation bin ^^ corresponds to the shell averaged correla- 
tion function 



?. 



^ ' C(OdV, 



V, 



(6) 



Vi 



where Vi is the volume of the shell. The covariance of this 
estimate is given by 



Akk'^-. 



dV / dVQ(r,r') 



27r2 



Joik,i)jo{k,j)ap{k), 



where 

jo{k,i) = TT I jo{kr)A^r. 



V, 



(7) 

(8) 

(9) 



Vi 



The dotted lines in each panel of Fig.|S]show the estimates of 
the variance for the various samples computed using Eq. (|8]). 
For the halo samples the variance is rescaled by the same 
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Figure 8. Comparison of the variance measured from our ensem- 
ble of simulations (solid line histogram) for the halo correlation 
functions in halo samples 1, 2 and 3 (lower panel, middle panel 
and upper panel respectively) with the predictions from Eqs. ^ 
(dashed lines) and llOII (doted lines). 



bias factors as the measured correlation functions. In all 
cases the variances measured from the ensemble and the 
theoretical predictions are in very good agreement. 

Fig. [7] shows a comparison of the correlation matrix 
Cor^{i,j) = C^{i,j)/a^{i)a^(j) predicted by Eq. (|8} with 
the estimate from the ensemble of simulations. The above 
diagonal parts of the matrices show the measured correla- 
tion matrices for the dark matter (panel (a)), and the haloes 
in samples 1, 2 and 3 (panels (b), (c) and (d) respectively). 
The lower triangular parts show the theoretical predictions. 
In general, Eq. (|HJ accurately describes the elements close to 
the diagonal, particularly for the dark matter. Further away 
from the diagonal, Cor(i,j) is noisy and the agreement is 
not so clear. 

ISmith et al.l (|2008r ) performed a similar comparison and 
reach the same overall conclusion as we do. These authors 
applied Eq. ((Sjl to describe the covariance matrix of the bin 
averaged correlation function in the limit Ar/r <C 1. In or- 
der to include the effect of the bin width, they re-wrote the 
contribution to the covariance coming from the term involv- 
ing 1/n^ as 



Vn2 






■jo{kri)ji)(krj) = 5i^. 



Vn^V,' 



(10) 



which corresponds to the exact result from Eq. (|8]) in the 
Poisson clustering limmit {P{k) <C 1/n)- This correction 
includes the effect of binning only in the last term of Eq. ([5} 
and ignores the effect of the shell-average in the first two 
terms. Fig. [8] shows a comparison of the variance measured 
from our ensemble of simulations (solid line histogram) for 
the halo correlation functions in halo samples 1, 2 and 3 
with the predictions from Eqs. ([8} and (|10|l . It can be clearly 
seen that even for the value of Ar/r = 0.03, Eq. (|10|) slightly 
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Figure 9. Upper panel: the linear perturbation theory power 
spectrum, Piin(fc) (dot-dashed line) compared with the "de- 
wiggled" power spectrum, Pdw(^) (solid line), as defined by 
Eq. I llll l. The individual terms on the right hand side of Eq. Illll l 
which make up the dewiggled power spectrum are also shown by 
short and long dashed lines as labelled. Lower panel: the ratios of 
the linear (dot-dashed line) and dewiggled power spectra (solid 
line) to a smooth, cold dark matter only power spectrumPnw(fc)- 



overpredicts the true variance, while the full expression from 
Eq. ^ gives a more accurate description. 

Eq. (|8]) is very useful for predicting the error and full co- 
variance matrix of two-point correlation functions estimated 
from galaxy samples. This is a valuable tool to use in fore- 
casts of the likely constraints attainable on cosmological pa- 
rameters from present and future galaxy surveys. 



5 MODELLING THE FULL SHAPE OF i{R) 

5.1 Non-linear evolution 

The nonlinear evolution of density perturbations changes 
the shape of the power spectrum and the correlation 
function due to the coupling between different Fourier 
modes. Nu merical simulations hav e been used to model 
this effect jEftathiou et al.l [l983: iHamilton et all Il99ll : 



IPeacock fc Dodds 
combination of the 



19941: iSmith et all 120031 ). Based on a 
Hamilton et al.l (119911 ) scaling relations 
and the halo model, ISmith et al.l l|2003l) proposed a fitting 



function to describe the effects of non-linear evolution on the 
shape of the matter power spectrum. This formula is able 
to reproduce the non-linear power spectra for pure CDM 
models to an accuracy of around 7% and is the basis of the 

commo nly used halof it code^ 

iHuff et al.l (|2007h and ICrocce fc Scoccimarrd (|2008l ') 
showed that the behaviour of the non-linear power spectrum 
predicted by halof it fails to describe an important distor- 
tion to the baryonic acoustic oscillations seen in numerical 
simulations. Besides the change in the overall shape of P{k), 



non-linear evolution washes out the acoustic oscillations 
by erasing the higher harmonic peaks (IMeiksin et al.l 19991: 
Eisenstein et al. 2005; Springel ct al.ll2005l : lEisenstein et al.l 
2006b ; Angulo et al. 2005 , 2008 ). This effect can be mod- 
elled by computing a damped or 'dewiggled' p ower spectrum 
iTegmark et all l2006l : lEisenstein et aLil2006bl ) : 



Pdw(fc) = Plin(fc)G(fc) + Pnw(fc)(l - G{k)), 



(11) 



where Piin(fc) is the linear theory power spectrum and 
Pnw{k) is a smooth, linear theory, cold dark matter only 
power spectrum, with the same shape as Piin(fc) but with- 
out any baryonic oscill ations (which can be com puted using 
the fitting formula of lEisenstein fc Hul l|l999l )'). The func- 
tion G(k) = exp [— (fc/\/2fcv,)^] regulates the transition from 
large scales {k <C A:*), where Pdw(fc) follows linear theory to 
small scales (fc ^ fc*) where the acoustic oscillations are 
completely damped. The upper panel of Fig. |9] compares 
Piin(fc) (dashed line) and Pdw(fc) (solid line) together with 
the different terms in Eq. (|lip . The damping of the BAG 
in the final power spectrum can be more clearly seen in the 
lower panel of Fig. [^ which shows the ratios of the linear 
and dewiggled power spectra to Pnw(fc). 

The degree of damping of the BAG is extremely sen- 
sitive to the value adopted for fc*. We treat fc* as a free 
parameter in the model for the form of the correlation func- 
tion. Some authors have attempted t o compute t he valu e 
of fc*. Based on the analysis of Eisenstein et al.l (|2006a ) , 
iTegmark et al.l l|2006r ) computed the value of the damping 
scale as a function of flm and the primordial amplitude of 
scalar fluctuations. As- These authors expressed their results 
in terms of ct* (where fc* — 1/ct*): 



soD{l + f) 



1/3 



As 



0.6841 



(12) 



Here, so = 12.4 h~^ Mpc is a reference scale, D is the growth 
factor (which is where f2m enters) , / = d log - D/dloga, and 
As fo llows the normalisation convention of ISpergel et al.l 

Eq. ITT} gives a purely phenomenological description of 
the damping of the acoustic oscillations found in numerical 
simulations. Using renormalized per turbation theory (RPT), 
ICrocce fc Scoccimarrd (|200q . |2008| ) give a theoretical justi- 
fication of the ability of Eq. (Illjl to describe non-linear evo- 
lution. According to RPT, the first term on the right hand 
side of Eq. (|lip describes the growth of a single mode, quan- 
tified by the propagator function G{k). In the high-fc limit 
the propagator is given by the Gaussian form with fc* given 
by (|Crocce fc Scoccimarro|[200d : lMatsubara II2008I ) 



fc* = 



^JdkPunik) 



-1/2 



(13) 



In the next section we shall compare the value of fc* obtained 
from Eqs. (I12|) and (|13|l with the best fitting value obtained 
by requiring our model for the correlation function to repro- 
duce measurements made from the N-body simulations. 

The second term on the right hand side of Eq. (|lip can 
be interpreted as the power generated by the coupling of 
Fourier modes on small scales, Pinc(fc). The term Pmc(fc) is 
negligible on large scales (small fc), but dominates the total 
power on small scales (high fc). For the scales relevant to 
the analysis of BAG, which are of the order of fc*, Pmc has 
a similar amplitude to Pnw(fc)(l — G(fc)). 
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Figure 10. An illustration of how the appearance of the correla- 
tion function changes after applying differing degrees of damping 
to the acoustic oscillations. The undamped, linear theory correla- 
tion function is shown by the solid line. The damped or "dewig- 
gled" correlation functions are shown by the dot-dashed lines, for 
different values of fc*. For small damping (large fc^), 5dw(^) de- 
viates by a small amount from 5iin('')- As the damping of the 
acoustic oscillations becomes stronger (smaller A;*), the peak is 
gradually erased and approaches the zero-baryon or "no-wiggle" 
correlation function, 5nw('") (dashed line) as fc* — > 0. 



Figure 11. Comparison of the z = real-space dark matter two- 
point correlation function averaged over the ensemble of simula- 
tions (open points) with: (i) the linear theory correlation function 
^lin(r) (solid line), (ii) an estimate of the non-linear correlation 
function §"'(?') computed using halof it without damping of the 
acoustic oscillations (dashed line), (iii) the dewiggled linear the- 
ory correlation function 5dw(^) defined by Eq. H14| l (dot-dashed 
line) and (iv) a dewiggled correlation function after being non- 
linearised using halof it C5j,,(»') (long-dashed line). The errorbars 
indicate the variance between the different realizations in the en- 
semble of simulations. 



Note, however, that according to RPT the propagator 
G{k) only behaves as a Gaussian in the high-fc limit. Be- 
sid es Pnw(k)(l — G(k)) is a smoo th function while, as shown 
by ICrocce fc Scoccimarrd l|2008t ). the term Pmc does show 
acoustic oscillations, although of a much smaller amplitude 
than P{k). This differences constitute the limitation of a 
modelling based on Eq. (|ll|l . We will return to this point in 

§E3] 

Taking the Fourier transform of Eq. ([TTJ , the two-point 
correlation function will be given by 



Cdw(r) = Ciin(r) ® G{r) + ^nw(r) ® (1 - Gir)), 



(14) 



where the symbol (g) denotes a convolution and G(r) is the 
Fourier transform of G{k). The information about the acous- 
tic oscillations is contained in the first term, which repre- 
sents the convolution of the linear theory correlation func- 
tion with a Gaussian kernel. This convolution implies that in 
the correlation function, the damping of the higher harmonic 
oscillations cause s the acou stic peak to broaden and shift 
to sm aller scales (ISmith et al. 2008; Croccc & Scoccimarro 
120081 ). Fig. [TD] shows the dewiggled correlation functions, 
^dw(^), obtained for different values of fc^ (dot-dashed lines). 
For weak damping (large fc*), only the highest harmonic 
oscillations are damped and ^dw('') deviates only slightly 
from ^iin(»') (solid line). As the damping becomes stronger 
(smaller fc*), the impact on the shape of the peak increases, 
with ^dvtir) approaching ^nw('') (dashed line) as fe* — > 0. 

A full description of the non-linear two-point correla- 
tion function must include both the overall change in shape 



of the correlation function and the damping of the acous- 
tic oscillations. In order to assess which of these effects is 
the more important for the analysis of the acoustic peak in 
^(r). Fig. 111! compares the mean z — real-space correla- 
tion function of the dark matter measured from the ensem- 
ble of simulations (open points) with the following models 
for the correlation function: (i) the linear theory correla- 
tion function 5iin(7') (solid line), (ii) a non-linear correla- 
tion function 5"'(r") computed using halof it, without any 
damping of the acoustic oscillations (short-dashed line) , (iii) 
the dewiggled linear theory correlation function £,dw{r), com- 
puted as described by Eq. (I14p (long-dashed line) and (iv) a 
dewiggled correlation function non- linearised using halof it 
Cdw('") (dot-dashed line). The errorbars indicate the vari- 
ance between the correlation functions measured from the 
different realizations in the simulation ensemble. 

It is clear from Fig.[TT]that the acoustic peak in the two- 
point correlation function at redshift z = shows strong de- 
viations from the predictions of linear theory. The halof it 
fitting formula fails to describe these deviations since its ap- 
plication only produces a small decrease in the amplitude of 
the peak without shifting its position. On the contrary, the 
linear theory dewiggled correlation function from Eq. (|14|l 
gives a very good description of the results of our numerical 
simulations, showing that the damping of the oscillations is 
the most important effect to include in the modelling of the 
real space correlation function on large scales. The incorpo- 
ration of the full change in shape of P{k) due to non-linear 
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Figure 12. Constraints on a and A;* obtained using the mean 
real-space (solid lines) and redshift-space (dashed lines) dark mat- 
ter correlation functions from the ensemble of simulations at red- 
shift 2: = 0, 0.5 and 1 (panels (a), (b) and (c) respectively). 



evolution produces very little difference in the shape of the 
acoustic peak in the correlation function. However, as we will 
see later (see section [5]3]) , this effect might be important on 
intermediate scales (r ~ 7Qh~^ Mpc). 
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Figure 13. The distribution of the values of a (upper panel) 
and fc* (lower panel) recovered by fitting the model of Eq. I I14II 
to the 2 = correlation functions measured from the different 
realizations in real (solid histogram) and redshift-space (shaded 
histogram). The dotted and dashed lines in the lower panel show 
the value of fc* predicted by Eq. II12I I and I I13I I respectively. 



5.2 The model in practice 

We now test whether or not the model for the correlation 
function described in § 15.11 returns unbiased constraints on 
cosmological parameters. In particular, we are interested in 
the constraints on the dark energy equation of state pa- 
rameter wde. We consider a very simple case in which we 
assume that the values of all cosmological parameters are 
known, apart from uidb, and analyse the constraints on this 
parameter. 

The change in the power spectrum due t o variations in 
wde in this simple case is described in detail bv lAngulo et al.l 
l|2008h . In order to measure the power spectrum of galaxy 
clustering, it is necessary to convert the angular positions 
and redshifts of the galaxies into comoving spatial separa- 
tions. This requires a choice to be made for the values of 
the cosmological parameters, including wq-e- The effect of a 
change in the value of wde from its true value to wde + ^iude 
is to modify the separations between pairs of galaxies, which 
leads to a change in the appearance of the power spec- 
trum. For small perturbations away from the true equation 
of state, the alteration in the measured power spectrum can 
be represented by a rescaling of the wavenumber from fctruc 
to fcapp. This change can be d e scribed by a 'stretc h' factor 
Q, defined by (JHuff et al.ll2"007l : lAngulo et al.ll2008l ) 



"^a 



kt^ 



(15) 



In the two-point correlation function, there will be an equiv- 
alent shift from scale rtruc to rapp — rtrue/a. We note that 
the variant models which arise from perturbing the equation 
of state in this way do not match observations such as the 
location of the Doppler peaks in the cosmic microwave back- 
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Figure 14. Panel (a): dewiggled correlation functions corresponding to A:* = 0.09, 0.1, 0.11, 0.15. 0.17 and 0.21, for a fixed stretch 
parameter, q = 1 (dot-dashed lines). To increase the dynamic range, we show the residuals g(r) — §dw('')i where ^dwC*") is the dewiggled 
correlation function for the best fit value of fc* = 0.115 /iMpc~^. The most important deviations between ^^^ and the measured 
correlation function is found at scales smaller than the position of the peak, with r < 90 h~^WLpc. Panel (b): Difference ^(r) — ^dwC*") for 
a = 0.978, 0.986, 0.996, 1.006, 1.014 and 1.022, fixing k^, = 0.115 hMpc"! (dot-dashed lines). A value of a slightly less than 1 re-scales 
the position of the peak in the model correlation function and improves the agreement with the measurements in small scales, causing 
the slight bias of the model thowards a < 1. The arrows indicate the effect of increasing fc<^ and a on the model correlation functions. 



The errorbars show the variance in the mean correlation function, 



that is (Tj : 
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ground (see lAngulo et al.ll2008n . Nevertheless, our primary 
goal here is to compare the constraints on the stretch param- 
eter obtained from the correlation function with those which 
result from the power spectrum, under the same idealized 
conditions. 

We analyse the constraints on the stretch parameter 
using the two-point correlation functions in real-space mea- 
sured from the different L-BASIIC II realizations, as de- 
scribed in Section 14.21 as well as from the mean correlation 
function from the ensemble. In doing so, we use only in- 
formation from the shape of the correlation function, and 
not its overall amplitude, which for real observational data 
can be affected by bias and redshift-space effects (see Sec- 
tions 15.41 and 15. 5p . We also consider fc* as a free param- 
eter, instead of fixing its value according to Eq. (|12|l . To 
explore this simple parameter space we constructed Monte 
Carlo Markov Chains (MCMC). The constraints obtained 
in this way for redshifts z — 0, 0.5 and 1 are summarised 
in Table 15.31 The range of pair separations used in the fit is 
60 < (r//i~^Mpc) < 180; the results are not sensitive in de- 
tail to the choice of limits in r. Unless otherwise stated, the 
quoted allowed ranges for the constrained parameters corre- 
spond to the 68% confidence level according to the variance 
from our ensemble of simulations. 

The solid lines in Fig. [l2]show the two-dimensional con- 
straints on a and fc* obtained using the mean dark matter 
correlation function in real-space measured from the ensem- 



ble of simulations. The constraints on fc* can be compared 
with the predictions of Eq. (|12p which for this cosmology 
gives fc* = 0.095 /iMpc~^ at z = 0. Here we find that at 
z = 0, the shape of the mean correlation function can be 
more accurately described by fc* — 0.115 ± 0.009 /iMpc~^, 
which is 2 — a away from the prediction of Eq. (|12p . On the 
other hand, this result is in excellent agreement with the 
prediction from Eq. (|13p which gives fc* = 0.117 ft Mpc~^ at 
z = 0. The lower panel of Fig. [TH] shows the histogram of 
the values of fc* obtained for the z — correlation functions 
of each L-BASIIC II realization, which is completely consis- 
tent with the constraints obtained from the mean correlation 
function of the ensemble. The dotted line shows the value 
of fc* computed using Eq. (|12|l which, as we mentioned be- 
fore, predicts a stronger damping of the acoustic oscillations 
than we find in the numerical simulations for this cosmol- 
ogy. This implies that if one applies Eq. (|12p , this may result 
in systematic errors in the values of the derived cosmologi- 
cal parameters. The dashed line shows the prediction from 
Eq. (|13p which is remarkably consistent with the mean value 
of fc* recovered from our ensemble of simulations. 



At higher redshift, the acoustic oscillations are less 
damped and higher values of fc* (corresponding to weaker 
damping) are favoured, with fe* = 0.140±0.010 feMpc"^ re- 
turned at 2 = 0.5 and fc* = 0.159 ± 0.013 ft Mpc"^ at z = 1. 
The constraints on fc* also broaden since, as can be seen in 
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Fig. |4l for higher values of fe* the shape of the acoustic peak 
becomes less sensitive to variations in this parameter. 

At redshift z = 0, the mean correlation function gives 
tight constraints on the stretch factor, with a = 0.996 ± 
0.006, showing that it is better described with a value of a 
slightly less than 1. The upper panel of Fig. 1131 shows the 
histogram of the values of a recovered from the different re- 
alizations, which is centred on the value a = 1, but which 
also shows a small tendency towards a < 1 which can be de- 
tected by the slight skewness of the distribution. At higher 
redshifts the position of the acoustic peak can be more pre- 
cisely determined since it is less affected by the damping 
of the higher harmonic oscillations. This results in tighter 
constraints on a, with a — 0.998 ± 0.004 at z = 0.5 and 
= 0.997 ±0.003 at 2 = 1. 

The origin of this slight bias towards a < 1 can be un- 
derstood by analysing the effect of variations in fc* and a on 
the model correlation function. Fig. I14f a') shows the corre- 
lation functions corresponding to different values of fc* for 
a fixed stretch parameter, a — 1 (dot-dashed lines). To in- 
crease the dynamic range, we plot the residuals ^{r) — ^dv{r), 
where ^dw('") is the dewiggled correlation function for the 
best fit value of fc* — 0.115 /i Mpc~^ . For this reference value, 
^dw accurately describes the shape of the acoustic peak in 
the mean correlation function measured from the ensemble 
of simulations. There are, however, differences between the 
model of Eq. (|14p and the measured correlation function 
both on scales smaller and larger than the acoustic peak 
which present a systematic structure. 

At scales around r ~ 150/i~^Mpc, Eq. ((TUl slightly 
overestimates the true value of the correlation function, 
which follows more closely the result from linear pertur- 
bation theory. On these scales, the value of ^dw is almost 
completely determined by the first term of Eq. (|14[) . with 
the second term giving only a slight correction. This differ- 
ence may indicate that the function G{k) that controls the 
damping of the oscillations is not exactly Gaussian, thereby 
changing the shape of the first term of Eq. (|14|) . or that the 
second term in this equation is not a good approximation 
to ^uic{r), and underestimates the full contribution of this 
term. 

The most important deviations between the model 
and the measured correlation function are found on scales 
smaller than the position of the peak, with r < 90 ft^^Mpc. 
On these scales the model underestimates the amplitude of 
the measured correlation function. 

Fig. 1141' b) shows the correlation functions obtained for 
different values of a, fixing fc* = 0.115 hMpc~^ (dot-dashed 
lines). As stated before, for a — 1 the model fails to repro- 
duce the measured correlation function. Using a value of a 
slightly less than 1 re-scales the position of the peak and 
increase the amplitude at small scales in the model correla- 
tion function in a way that gives better agreement with the 
measurements from the simulations, causing the slight bias 
of t he model towards a < 1 . 

ICrocce fc Scoccimarrd l|2008l ) showed that (,mc{r) intro- 
duces an additional shift in the acoustic peak which is not 
taken into account by the second term in the right hand side 
of Eq. (|14|l . At z = this shift can lead to a bias towards 
Q > 1 of the order of 0.8% if this parameter is fitted using 
only the position of the peak (see their Figure 7). This ap- 
parent contradiction with our result of a small bias towards 



a < 1 is due to the fact that our constraints are based on the 
full shape of the correlation function on large scales, and not 
only on the position of the peak. In this case, as we shown 
above, the biggest difference between the models and the 
measured correlation functions is not given at the scales of 
the peak, but at smaller scales, where non-linear evolution 
increases the amplitude of ^(r) with respect to the predic- 
tion of Ea. (|14p . This difference is responsible for the shift in 
the preferred value of a. 

5.3 Improving the model 

Although the model described by Eqs. ((TT)) and H14p gives a 
good description of the distortions in the shape of the peak 
in the correlation function due to the damping of the acous- 
tic oscillations, it is slightly biased towards a < 1. This bias 
is smaller than the sample variance predicted in the corre- 
lation function for the volume of the L-BASIIC simulation 
or expected in ongoing galaxy surveys. However, the much 
smaller sample variance anticipated in future galaxy surveys 
such as Euclid (|Robberto et al. IbOOTl ) or ADEPT, will mean 
that this level of inaccuracy in the modelling of the corre- 
lation function may lead to the introduction of systematic 
errors in the recovered dark energy equation of state param- 
eter u;de. In this section we explore different ways to correct 
for this small bias. 

The model of Eq. (|ll|l can be extended to account for 
the change in the overall shape of the power spectrum due 
to non-linear evolution by 



PSUk) 



1 + ' 



1 + Ak + Bk'- 



PiUk) = f{k)Pi„{k). 



(16) 



The factor /(fc) could also be used to model a scale de- 
pendent bias factor. Thi s model for non-l inear evolution is 
based on the Q-model of ICole et al.l l|2005h . modified by the 
addition of a new parameter, B, with the aim of achieving a 
better description of the behaviour of the non-linear power 
spectrum at high k. The value of the new parameter B de- 
pends on the choice of Q. In practice, we fixed the value of B 
using the relation B = Q/10, which gives the approximate 
behaviour of the non-linear power spectrum at large fc, al- 
though the correlation function is almost insensitive to this 
choice. This scheme is an alternative means to describe the 
non- linear growth to using halof it. This expression can be 
Fourier transformed to obtain a non-linear correlation func- 
tion ^dw('')i whose shape is only weakly dependent on the 
precise value of B. 

Fig. I15f a) shows the quantity ^JJw — (,\® G for differ- 
ent values of the parameter Q fixing a = 1 and ki, = 
0.115 ft Mpc"^ (dot-dashed lines). For Q = 13 (thick dot- 
dashed line), the non-linear correction factor from Eq. p6|) 
changes the shape of the model correlation function, improv- 
ing the agreement with the measurements from the simula- 
tions on scales r < 90 ft~^Mpc. On large scales the measure- 
ments follow more closely the predictions from linear per- 
turbation theory and the models overpredict the amplitude 
of the correlation function. 

Table 15.31 also shows the constraints on a and fc* ob- 
tained by including the scale dependent correction factor 
defined by Eq. (|16p , on allowing Q to vary as a free parame- 
ter. As the model correlation function is fairly insensitive to 
the value of Q the constraints on this parameter are weak. 
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Figure 15. Panel (a): The quantity (^jj^ - ^iin G) for Q = 1, 7, 13, 19 and 26, fixing a = 1 and fc* = 0.115/iMpc-i (dot-dashed lines). 
For Q = 13 (thick dot-dashed line), the non-linear correction factor from Eq. JTBJ changes the shape of the model correlation function, 
improving the agreement with the measurements from the simulations on scales r < 90 h~^y[pc. Panel (b): The quantity (^"' — ^n^ (g) G) 
computed using Eq. ^9^ with Amc = 0.15, 0.25, 0.37, 0.45 and 0.55, fixing « = 1 and fc* = 0.115 /iMpc-^ (dot-dashed lines). This 
modelling gives a better description of the residuals of the measured correlation functions with respect to (^iin ^ G) close to the acoustic 
peak. The solid lines on both panels show the difference (^Jw ~ ?lin ^ G), that is, the second term on the right hand side of Eq. II14I I. 
The arrows indicate the effect of increasing Q and Amc on the model correlation functions. The errorbars show the variance in the mean 
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However, there is a tendency towards lower values of Q with 
increasing redshift. As can be seen in Fig. 1151 varying Q 
changes the amplitude of the acoustic peak. This is why 
when this parameter is included in the analysis, the best fits 
shift towards lower values of fe*. This reflects the fact that 
when the non-linear correction of Eq. (|ll|l is not included, 
the value of fc* recovered from the correlation function con- 
tains some information on the distortion to the overall shape 
of -P(fc) and not just the way in which the oscillations are 
damped. 

The constraints on a are essentially unchanged upon the 
incorporation of the non- linear shape correction. The best fit 
results, as before, show a small bias towards q < 1. Hence, 
although the model presented in Eq. (I16|) gives the freedom 
to correct for differing degrees of non-linear evolution or even 
a scale dependence of bias, it is still unable to describe the 
form of the measured correlation functions to correct for this 
small systematic shift in a. 

The solid curve in Fig. \TS\ shows the quantity ^2w ^ 
^1 ® G, that is, the second term on the right-hand side of 
Eq. ((14} . As pointed out in § [S] this term is a smooth func- 
tion with n o information about th e acoustic peak. On the 
other hand, ICrocce fc Scoccimarrd (J2008) showed that the 
term ^mc(r) does contains information about the acoustic 
scale. This fact imposes a limitation on the ability of a model 
based on Eq. (|14|) to reproduce the full shape of the non- 
linear correlation function. 



Using stand a rd perturbation theory, 

ICrocce fc Scoccimarrd (|200a ) showed that on the scale of 
the acoustic peak the main contribution to ^mc('') will be 
of the form 



Cmc(r)(xe,'i„e,\i'(r 



(17) 



where ^ii„ is the derivative of the linear theory correlation 
function and 



^l^(r) 



v-ieii„( 



/ 



47r / Piin(fc) ji(fcr)fcdfc. 



(18) 



Based on this result, Crocce & Scoccimarro proposed an 
approximate phenomenological model to describe the effect 
of the mode coupling on the correlation function near to 
the acoustic scale, in which ^mc(»') ~ ^mcCun (since Cun ('') 
varies little on these scales). We tested this ansatz by com- 
paring the results from our simulations against the non- 
linear correlation function given by 



Cni(r) = Cim(r) ® G(r) + ^„ 



= CC^(r), 



(19) 



where y4,nc is a free parameter. Fig. I15f b) shows the non- 
linear correlation functions obtained for different values of 
Amc fixing a = 1 and fc* = 0.115 /iMpc~^ (dot-dashed 
lines). It can be seen that this model can describe the shape 
of the residuals of the measured correlation functions with 
respect to ^un ® G close to the acoustic peak. At smaller 
scales, where the approximation is not so accurate, the 
model underestimates the correlation functions. 
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Table 2 also shows the results obtained by applying 
this modelling. The constraints obtained on fe* are slightly 
smaller than those obtained using Eq. (|14|l . The most impor- 
tant efTect of the incorporation of the second term of Eq.([l9]) 
is that it helps to alleviate the problem of the small bias to- 
wards a < 1 found in § [S^] yielding a = 1.003 ± 0.008 at 
z — 0, with a slight increase in the allowed region. At higher 
redshifts the constraints are tighter, with a — 1.002 ± 0.005 
a,t z — 0.5 and a = 1.000 ± 0.003 at z = 1, in better agree- 
ment with a = 1. The ability of this simple ansatz to im- 
prove the obtained constraints might indicate that the im- 
plementation of a full calculation of ^mc using RPT over the 
full range of scales included in the analysis could help to 
improve the constraints even further. 



5.4 Redshift-space effects 

When dealing with spectroscopic surveys, the distance to a 
galaxy is inferred from its measured redshift. As the red- 
shift is affected by the peculiar motion of the galaxy along 
the line of sight, the radial distance so obtained does not cor- 
respond to the true distance to the galaxy. This changes the 
clustering pattern of galaxies, leading to differences between 
the power spectrum and correlation function measured in 
redshift-space, Ps{k) and ^s{r), and their true real-space 
counterparts. These differences are called redshift-space dis- 
tortions. As the peculiar velocity field of galaxies is induced 
by the underlying matter distribution, it is possible to model 
the distortions that will arise in a given cosmological model. 
On large scales where linear theory is applicable, the 
peculiar velocities take on the form of coherent bulk flows 
towards overdense regions. This leads to an increase in the 
amplitude of Ps{k) and ^s{r) compared to the real-space 
values. Under the assumption of linear perturbation theory 
and the plane parallel approximation, the scale-independent 
boost in power produced by this effect is given by (jKaiserj 
1 19871 ) 



P{k) 



^S 



fi-hfcV) 



(21) 
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C(r) 



5'^+ 5^' 



(20) 



where /3 = f/b and b is the bias factor (which for the dark 
matter is simply 6=1). The same factor applies for the 
power spectrum. 

On small scales peculiar velocities are dominated by 
the random motions inside virialized structures. This makes 
bound structures such as dark matter haloes to elongated 
when mapped in redshift-space, an effect commonly known 
£18 'fingers of god'. This smearing of structure causes a damp- 
ing of the power spectrum and the correlation function on 
small scales. A complete description of redshift space distor- 
tions must include both of those regimes (for an empirical 
des cription, sec P eacock & Dodds , 1994). 

IScoccimarro I (|2004l ') predicted that in the case of the 
power spectrum, the boost factor of Eq. (|20|) is reached 
asymptotically on very large scales {r > 100/i~^Mpc) and 
confirmed this with intermediate scale N-body simulations. 
Similar conclusions were reached with much larger simula- 
tions by Angulo et al. (2008). For smaller scales the power 
is damped by r andom motions in a way that can be ap- 
proximated by (jPeacock fc Dodds I Il994l : iPark et all 1 1994 
ICole et al.|[l994l ) 



where cr is a free parameter connected to the pairwise veloc- 
ity dispersion. This implies that in order to use information 
contained in the full shape of the power spectrum to obtain 
constraints on cosmological parameters, a detailed model of 
redshift space distortions must be implemented. 

The open triangles of Fig[SJa) show the mean two point 
correlation function measured in redshift-space from the en- 
semble of simulations, rescaled by the boost factor S from 
Eq. (|20|l . This accurately corrects the amplitude of ^s (r) to 
match that of the real-space ^(r) over the range of scales 
plotted. The most important effect on the correlation func- 
tion of the decrease in power on smaller scales (as described 
by the extra factor in Eq. (|21[) ') is to produce an extra damp- 
ing of the acoustic oscillations. This can be accurately mod- 
elled by a smaller value of fe*. 

Hence the treatment of redshift-space distortions in the 
correlation function is somewhat simpler than for the power 
spectrum, at least for the relatively large scales we consider 
here. If the amplitude of the correlation function is marginal- 
ized over when constraining cosmological parameters, the 
presence of the boost factor S is irrelevant. In this case, 
redshift-space distortions simply produce a stronger damp- 
ing of the acoustic oscillations, which is reflected in the flt 
to the measurements preferring smaller values of fc*. This 
can be seen in Fig. 1121 which shows the two-dimensional 
constraints on a and fc*. In redshift-space, the measure- 
ments from the simulation prefer lower values of fc* than in 
real-space. The constraint on fc* is also tighter in redshift- 
space, with k^, = 0.097 ± 0.004 ft Mpc"^ at 2 = 0. This can 
also be seen in the lower panel of Fig. 1131 which shows the 
distribution of the values of fc* recovered from the z = 
redshift-space correlation function of the different realiza- 
tions (dashed histogram). The best value for fc* in redshift- 
space is in better agreement with the prediction of Eq. (|12p . 
The same situation is found at higher redshifts, where the 
preferred values for fc* are fc* = 0.107 ± 0.005 ft Mpc~^ at 
z = 0.5 and fe* = 0.124 ± 0.005 ft Mpc"^ at 2 = 1, in both 
cases lower than the corresponding values in real-space. 

The constraints on the stretch factor found using 
redshift-space clustering are weaker than in real-space, with 
a = 1.002 ± 0.009 a,t z — 0, showing a considerable in- 
crease in the width of the allowed region for this parameter. 
However, this result does not show the small bias seen in 
a in real-space, since redshift-space distortions alleviate the 
problem of non-linear evolution at small scales. At higher 
redshift, the preferred values again towards a < 1, with 
a = 0.995±0.006 at z = 0.5 and a = 0.994±0.005 a.t z = 1, 
with a = 1 on the edge of the 68% c.l. regions. 

When the non-linear correction of Eq. (|16|l is applied, 
the constraint on a remains unaltered. As when using real- 
space measurements, the same trend towards lower values 
of Q with increasing redshift is found, and the mean values 
for fc* are lower than the ones found without applying this 
non-linear correction. 

Table 2 also shows the constraints obtained comparing 
the model of Eq. (I19|) against the redshift-space dark mat- 
ter correlation functions. The mean values for fc* recovered 
in this case are in excellent agreement with those obtained 
applying Eq. (|14|l . On the other hand, the constraint on a 
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Sample id 



fe2 



Real-space 
ki,/{hMpc-^) a 



Redshift-space 
ki,/{hMpc~'^) a 





0.0 


1.00 


1.34 


0.115 ±0.009 


0.996 ± 0.006 


0.097 ±0.004 


1.002 ± 0.009 


DM 


0.5 


1.00 


1.56 


0.140 ±0.010 


0.998 ±0.004 


0.107 ±0.005 


0.995 ± 0.006 




1.0 


1.00 


1.69 


0.159 ±0.013 


0.997 ±0.003 


0.124 ±0.005 


0.994 ±0.005 




0.0 


2.81 


1.19 


0.117 ±0.017 


0.997 ±0.012 


0.094 ±0.013 


0.992 ±0.012 


Halos 1 


0.5 


5.57 


1.21 


0.146 ±0.024 


0.995 ±0.011 


0.107 ±0.014 


1.001 ±0.013 




1.0 


11.58 


1.17 


0.22 ±0.10 


0.997 ±0.014 


0.131 ±0.028 


0.997 ±0.017 




0.0 


3.40 


1.17 


0.170 ±0.065 


0.997 ±0.010 


0.099 ± 0.024 


0.998 ±0.020 


Halos 2 


0.5 


7.35 


1.18 


0.20 ±0.10 


1.006 ±0.015 


0.160 ±0.066 


1.002 ±0.020 




1.0 


16.09 


1.15 


0.27 ±0.15 


1.003 ±0.022 


0.22 ±0.13 


0.998 ±0.019 




0.0 


6.46 


1.12 


0.121 ±0.023 


0.994 ±0.012 


0.109 ±0.019 


0.990 ±0.015 


Halos 3 


0.5 


13.14 


1.13 


0.203 ±0.076 


0.988 ±0.017 


0.187 ±0.070 


0.983 ±0.017 




1.0 


28.13 


1.11 


0.28 ±0.13 


1.000 ±0.025 


0.24 ±0.13 


1.000 ±0.021 



Table 1. Results of applying the general fitting procedure described in Sec. 15.11 to the correlation function of the different samples 
analysed in real-space and redshift-space. The first column gives the label of the sample, as defined in Section 2. The second column 
gives the redshift output. Column 3 gives the real-space bias factors computed as described in Sec. 15.51 Column 4 gives the value of the 
redshift-space boost factor computed according to Eq. I|20ll . Columns 5 and 6 give the constraints on ki, and a obtained by fitting ^dw 
as in Eq. I I14I I in real-space while columns 7 and 8 show the results obtained using the redshift-space measurements. All quoted values 
correspond to the mean and 68% c.l. of the respective parameter. 



show much better agreement with the value a = 1, espe- 
ciaUy at higher redshif t s. Thi s implies that, as suggested by 
ICrocce fc Scoccimarrd (2008), the ansatz given by Eq. (|19p 
can also correct for the small bias towards a < 1 when using 
redshift-space information. 



5.5 Halo bias 

According to the current picture of galaxy formation, galax- 
ies are hosted within dark matter halos. Hence, a model 
for the clustering of dark mater halos is a necessary first 
step toward s an understan ding of galaxy bias. Moreover, as 
discused in lAngulo et al.l (|2005i 'l , galaxy clusters can be a 
useful tracer of BAO in their own right, offering an impor- 
tant cross-check of the results obtained from galaxy samples. 
Clusters are more strongly clustered than galaxies, thereby 
increasing the amplitude of the correlation function signal 
on the very scales relevant to BAO, on which the ampli- 
tude of the typical galaxy correlation function is very low. 
Future cluster catalogues will map in an homogeneous way 
volumes large enough for an accurate determination of the 
signa l of the acoustic peak IjRuhl et al. ll2004l : lBartlett et al I 
l200g ). For these reasons the analysis of the way in which the 
imprints of the acoustic oscillations on the halo two-point 
correlation function change as clustering evolves is of great 

imp ortance. 

ISmith et al.l (|2008l ) performed a detailed analysis of the 
halo correlation function on the scales relevant to studies of 
acoustic oscillations. They focused on the determination of 
the existence of shifts in the peak of the correlation func- 
tion from the predictions of linear perturbation theory. If 
the peak in the halo correlation function is identified with 
the sound horizon at recombination in order to constraint 
the dark energy equation of state parameter, these shifts 
introduce important biases on the obtained constraints. As 



we have discused in sec. [5] the position of the peak in the 
correlation function should not be considered as a precise 
cosmological tool since it does not correspond to the ex- 
act value of the sound horizon s. Instead, a detailed mod- 
elling of the f ull sh ape of the correlation function is required. 
ISmith et al.l l|2008h also proposed an analytic model for the 
halo two point correlation function in real-space based on so- 
lutions to the pair conservation equation using characteristic 
curves and a model of the pairwise velocity distribution of 
the haloes. This model is able to reproduce the main trends 
found in numerical simulations, showing that it contains the 
most relevant physical processes that control the deviations 
from linear theory. Here we follow a more practical approach 
and analyse the biases or systematic effects introduced in the 
constraints on cosmological parameters if the modelling of 
sec. 15.11 is applied to describe the shape of the correlation 
function of different halo samples. 

On the very large scales, where halo biasing is determin- 
istic, a simple linear relation is expected to hold between the 
correlation function of dark matter halos, S.hir), and the lin- 
ear theory matter correlation function, with an effective bias 
factor, 6^g = ^h/Ciin, independent of scale. This effective bias 
can be computed by taking the weighted average of the bias 
b{m) as a function of halo mass m, over the selected halo 
sample: 



/„°° ^p{m) h{ra)n{rn) Am 
J^ tp{Tn)n{m)dm 



(22) 



where n{m) is the mass function of dark matter halos, which 
gives the space density of halos in the mass interval m to 
+ 6m (e.g . iPress fc Schechter I Il974 iBond et al] Il99ll: 



Bower I [l99 
2001I1. b( 



Sheth. Mo fc Tormen 
ives the bias as a function of halo 



200 ll : IJenkins et al] 



Mo fc White Ill996l: IMo. Jing fc White Ill996l: Jjing Ill998l : 
ISheth. Mo fc Tormen I I2001I : ISeliak fc Warren llJOM ) and 
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Real-space 












Fits with £,f^ 






Fits with CSjc 




Sample id 


z 


K/{hMpc-^) 


a 


Q 


k^,/{hMpc''^) 


a 


^mc 




0.0 


0.107 ±0.013 


0.996 ± 0.006 


18 ±12 


0.1135 ±0.0089 


1.003 ±0.008 


0.37 ±0.22 


DM 


0.5 


0.126 ±0.018 


0.997 ±0.004 


19 ±12 


0.133 ±0.010 


1.002 ±0.005 


0.39 ±0.21 




1.0 


0.155 ±0.030 


0.997 ±0.003 


13 ±10 


0.155 ±0.013 


1.000 ±0.003 


0.27 ±0.18 





0.0 


0.115 ±0.021 


0.996 ±0.012 


16 ±11 


0.120 ±0.016 


1.002 ±0.013 


0.35 ±0.25 


Halos 1 


0.5 


0.151 ±0.032 


0.993 ±0.011 


14 ± 11 


0.142 ±0.022 


0.999 ±0.012 


0.41 ±0.28 




1.0 


0.20 ±0.10 


0.995 ±0.013 


15 ±12 


0.199 ±0.087 


1.003 ±0.015 


0.52 ±0.36 




0.0 


0.154 ±0.055 


0.997 ±0.011 


16 ±11 


0.158 ±0.050 


1.004 ±0.011 


0.63 ± 0.46 


Halos 2 


0.5 


0.20 ±0.10 


1.006 ±0.016 


15 ±10 


0.186 ±0.085 


1.010 ±0.016 


0.52 ±0.39 




1.0 


0.27 ±0.15 


1.004 ±0.021 


15 ±10 


0.27 ±0.16 


1.009 ±0.022 


0.58 ±0.42 



Halos 3 



0.0 


0.115 ±0.022 


0.994 ±0.012 


17 ± 12 


0.121 ±0.020 


1.003 ±0.015 


0.52 ±0.35 


0.5 


0.21 ±0.11 


0.988 ±0.017 


14 ± 11 


0.185 ±0.062 


0.997 ±0.018 


0.66 ±0.46 


1.0 


0.27 ±0.14 


1.000 ±0.025 


15 ±10 


0.28 ±0.15 


1.010 ±0.027 


0.78 ±0.59 



Redshift-space 



Fits with £f 



Sample id z ki,/{hh[pc ^) 



Q 



Fits with 5 



K/{hMpc-^) 



nl 
mc 



0.0 0.092 ±0.008 1.002 ±0.009 18 ± 12 

DM 0.5 0.102 ±0.010 0.994 ±0.006 16 ± 11 

1.0 0.122 ±0.013 0.993 ±0.005 13 ± 10 



Halos 3 



0.0969 ± 0.0051 
0.1074 ±0.0057 
0.1248 ±0.0056 



1.009 ±0.010 
1.003 ±0.008 
0.998 ± 0.006 



0.34 ±0.21 
0.29 ±0.19 
0.22 ±0.15 





0.0 


0.091 ±0.013 


0.991 ±0.012 


18 ±13 


0.096 ±0.012 


0.996 ±0.013 


0.37 ±0.26 


Halos 1 


0.5 


0.104 ±0.016 


1.000 ±0.012 


18 ±13 


0.110 ±0.012 


1.008 ±0.013 


0.52 ±0.33 




1.0 


0.127 ±0.028 


0.997 ±0.017 


16 ±12 


0.133 ±0.027 


1.004 ±0.018 


0.53 ±0.35 




0.0 


0.094 ±0.021 


0.996 ± 0.020 


17 ± 12 


0.099 ±0.021 


1.003 ±0.019 


0.53 ±0.39 


Halos 2 


0.5 


0.142 ±0.050 


1.001 ± 0.020 


16 ±12 


0.143 ±0.046 


1.008 ±0.020 


0.57 ±0.43 




1.0 


0.19 ±0.10 


1.000 ±0.020 


15 ±11 


0.20 ±0.12 


1.000 ±0.019 


0.50 ±0.37 



0.0 


0.104 ±0.019 


0.989 ±0.015 


18 ±13 


0.109 ±0.017 


1.000 ±0.016 


0.61 ±0.41 


0.5 


0.182 ±0.080 


0.982 ±0.017 


16 ±12 


0.150 ±0.045 


0.997 ±0.021 


0.93 ±0.62 


1.0 


0.23 ±0.14 


1.000 ±0.021 


15 ±11 


0.20 ±0.10 


1.016 ±0.028 


0.64 ±0.50 



Table 2. Results of applying the general fitting procedure described in Sec. 15.31 to the correlation function of the different samples 
analysed in real-space (top) and redshift-space (bottom). The first column gives the label of the sample, as defined in Section 2. The 
second column gives the redshift output. Columns 3, 4 and 5 show the constraints on fc*, o and Q obtained by fitting ^Jj' , with the 



non-linear shape correction from Eq. I|16| l. Columns 6, 7 and 8 show the constraints on fc* 
Eg. 1191 All quoted values correspond to the mean and 68% c.l. of the respective parameter. 



a and Amc obtained using the model of 



tpirn) represents the mass selection function applied to con- 
struct the halo sample. For the cases analysed here, %l){rn) is 
simply given by the limits of the mass bins, i.e. tp{m) = 1 if 
mi < m < m2 and tp{m) = otherwise. 

The effects of halo bias on the two-point correlation 
function can be seen in Fig. [S] The correlation functions of 
the different halo samples can be rescaled using a constant 
bias factor, b^, to agree remarkably well with the measure- 
ment for the dark matter. This means that the halo bias 
is independent of scale for the pair separations plotted. For 
each halo sample and redshift analysed, we computed the 
bias factors b^ that maximize the likelihood of the best fit 
^dw model to the dark matter correlation function. The re- 
sults are shown in Table 15.31 For z — 0, these values are 



6 — 1.68, 1.84 and 2.54 for the halo samples 1, 2 and 3 respec- 
tively. We used these bias values to rescale the halo-halo cor- 
relation functions in Fig. [S] and to compute the theoretical 
variances for each sample as described in Section [4.31 These 
values can be compared with those predicted by Eq. (I22JL 
Using the recipe for the mas s function from [Je nkin s ct^.l 
((2OO1I ) and the bias function of lSheth. Mo fc Tormeni (l200lh 
(with the modified parameters of lSheth fc Tormen II2OO2 I) we 
find that Eq. (|22|) gives a good description of the bias factors 
with 6cft — 1.62, 1.89 and 2.53 f or the same halo samples . 
Using the bias prescription from ISeliak fc Warren I (|200J), 
the dis agreement wit h the e stimates from the simulations is 
bigger. lAngulo et al.l l|2008h measured the bias factor as a 
function of halo mass directly from the BASICC simulation 
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Figure 16. Constraints on a and fc* obtained using the mean real-space (solid lines) and redshift-space (dashed lines) halo-halo 
correlation functions for the halo samples 1, 2 and 3 (first, middle and last column respectively) from the ensemble of simulations at 
redshift z = 0, 0.5 and 1. (upper, middle and lower rows respectively). 



and found good agreement with the analytic prescription we 
use here, for the modest peak heights corresponding to our 
halo samples. 

In order to analyse any possible biases or systematic 
effects that can be introduced when using information from 
halo samples to constrain cosmological parameters, we ob- 
tained constraints on fc* and a by applying the model de- 
scribed in Section [5] to our measurements of the halo two- 
point correlation function. Our results are summarised in 
Table 15.31 Besides the overall constant bias factors, which 
are ignored (marginalized over) by our analysis procedure, 
the estimates of S,hir) show a change in the shape and posi- 
tion of the acoustic peak with respect to ^iin(''), similar to 
the changes seen in the dark matter correlation function. 

Fig. 1161 shows the two-dimensional constraints on a and 
fc* obtained using the mean halo-halo correlation functions 
for halo samples 1, 2 and 3 in real (solid lines) and redshift- 



space (dashed lines) measured from the ensemble of simula- 
tions. As a consequence of the larger shot noise, the allowed 
regions for these parameters increases with respect to the 
ones obtained using the dark matter correlation function. 
The upper panel of Fig. ll7l shows the histogram of the values 
of a obtained for the z = halo-halo correlation functions 
for halo sample 1 of each L-BASIIC II realization, which 
is completely consistent with the constraints obtained from 
the mean correlation function of the ensemble. This distri- 
bution is consistent with q = 1 though it presents the same 
small tendency towards values of a < 1 found in § 15.21 with 
the dark matter. 



ISmith et al.l (J2008l ) performed a similar analysis, fitting 
the halo correlation function for different halo samples with 
the linear theory ^(r) smoothed with a Gaussian filter. They 
found a slight tendency towards larger shifts in the position 
of the acoustic peak, which implies smaller values of fc* , with 
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Figure 17. The distribution of the values of a (upper panel) and 
ki, (lower panel) recovered by fitting the model of Eq. 1141 1 to the 
2 = correlation function of halos in sample 1 measured from 
the different realizations in real (solid histogram) and redshift- 
space (shaded histogram). The dotted and dashed lines in the 
lower panel show the value of fc* predicted by Eq. II12I I and 1131 1 
respectively. 



increasing mass. We do not reproduce such a tendency here. 
The lower panel of Fig. [T7] shows the histogram of the val- 
ues of fc* recovered from each realization. As a reference, 
we also show the theoretical predictions from Eqs. ((12} and 
H13p. Although the peak of the obtained distribution for this 
parameter is consistent with the one found for the dark mat- 
ter, there is a tail towards higher values of fc*. The same 
behaviour is seen in the results obtained for the remaining 
halo samples. This tail causes the obtained values for fc* to 
be larger than for the dark matter, indicating that on av- 
erage the halo BAO experience weaker damping than we 
find for the dark matter. This tail is more important for 
the higher mass bins which then creates a tendency towards 
higher mean values of fc* with increasing halo mass, even 
though the peak of the likelihood functions for this param- 
eter are always consistent with the ones found for the dark 
matter. 

As we saw with the dark matter, the incorporation of 
Q as a free parameter of the model does not alter the con- 
straints on a. The values of Q preferred by the halo mea- 
surements decrease at higher redshifts, when non-linear dis- 
tortions are not as strong as at 2; = 0. On the other hand, 
the constraints on the stretch parameter obtained by apply- 
ing Eq. (|19|l are completely consistent with a — 1 showing 
that this modelling is able to correct the small bias towards 
a < 1 present in the previous results also in the case of halo 
samples. 

The amplitude of the halo correlation function, like 
the matter correlation function, is affected by the redshift 
space distortions. These effects can be correctly described 
by Eq. ((20} with the value of S computed using the bias 



Figure 18. Upper panel: A comparison of the real-space dark 
matter power spectrum averaged over the simulation ensemble 
(open points) with the linear theory power spectrum (dot-dashed 
line), the 'dewiggled' power spectrum from Eq. 1111 (solid line), 
and its non-linear version from Eq. II6I I (dashed line) computed 
with Q = 13. The dotted lines indicate the variance on P{k) esti- 
mated from the ensemble. Lower panel: The ratio of these power 
spectra to Pnw(fc). The results from the numerical simulations 
have also been divided by the non-linear growth factor given by 
Eq. ITet . 



factors described above for the different halo samples. The 
values of S obtained this way are shown in Table 2 and were 
used to rescale the redshift-space halo correlation functions 
plotted in Fig. [6| As we found for the dark matter, besides 
this increase in the amplitude, the halo correlation function 
in redshift-space shows a stronger damping of the acoustic 
signal than in real-space. This is reflected in the allowed re- 
gion for ki,, which shows a preference for lower values than 
are obtained for real-space data. The constraints on a show 
similar behaviour to those in real-space, with some deviation 
from Of = 1, which, again, is unaltered upon the inclusion 
of the non- linear correction of Eq. ((16} . As in the previous 
cases, this small bias is corrected by the implementation of 
the modeUing of Eq. ((T^ . 



6 POWER SPECTRUM VS CORRELATION 
FUNCTION 

The next generation of galaxy surveys will cover signif- 
icantly larger volume s than current surveys (e.g. Euclid, 
iRobberto et al. II2007I . or ADEPT). For the first time, sys- 
tematic effects on the scale of the BAO will become compa- 
rable to the sampling errors. The power spectrum and corre- 
lation function of clustering are affected in different ways by 
effects such as non- linear evolution, bias and redshift space 
distortions. In this section, we compare the constraints on 
the dark energy equation of state from these two statistics. 
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to assess which one yields the smallest random and system- 
atic errors. We compare the constraints obtained by fitting a 
model to the measured correlation function to those derived 
usi ng the power spectr um in the same simulations as used 
bv lAnguloet alj l|2008h . 

We first recap the power spectrum and correlation func- 
tion approaches to constraining the equation of state to high- 
light where differences may lie between them. 

Various power spectrum approaches have been proposed 
(e.g. Blake fc Glazebrook I 20031: Dolnev. Jain fc Takada 



200e 



2007 



Wang 



nr 



Eisenstein 



2006"; 'Porcival et al."2007a'; 'Soo 
Seo et a l. 2008). Angulo ct al. ( 2008) developed the 



method introduced by IPercival et al.l l|2007ah . in which no 
information is used from either the amplitude or large scale 
shape of the power spectrum. A reference spectrum is de- 
fined from the measured power spectrum, by applying a 
spline fit to the spectrum after rebinning into coarser bins. 
The measured power spectrum is then divided by this refer- 
ence spectrum to form a ratio which emphasises the appear- 
ance of the BAO. The philosophy behind this conservative 
approach is that dividing the measured power by a reference 
which is defined directly from the measurement reduces the 
impact of any long wavelength gradients in P{k) induced 
by non-linear effects or redshift-space distortions, at the ex- 
pense of losing some of the information contained within 
the power spectrum. Such gradients do exist in the bias fac- 
tor and redshift space distortions measured in simulation s 
l|Scoccimarro 112004 [sinith et al.l [20081 : 1 Angulo et all 120081 ). 
A further advantage of taking a ratio of power spectra is that 
this removes the need for an accurate model of such effects. 
The measured ratio is then compared to a ratio generated 
from a linear perturbation theory spectrum. In this case, 
a new reference spectrum is generated for each linear the- 
ory -P(fc). The linear theory ratio is 'de- wiggled' in a similar 
way to Eq. (Illf) . Angulo et al. used the wavenumber range 
fc = - 0.4/iMpc"^ in their fits. 

In the correlation function approach described in Sec- 
tion O we model the full shape of ^(r). The key step is the 
damping or "de-wiggling" of the BAO in the power spec- 
trum, as described by Eq. (|ll|l . before taking the Fourier 
transform of the power spectrum to obtain the correlation 
function. We also apply two empirical non-linear distortions 
to the form of the power spectrum, using Egs . (I16p and 
(1191) . The first approach, based on the Q-model of I Cole et al.l 
(^OOS ) , has little impact on the quality of the fit on the scales 
used, r = 60-180 h ~^ Mpc. On the other hand t he second 
model, proposed bv lCrocce fc Scoccimarrd l|2008r ) based on 
an ansatz inspired by RPT, gives improved constraints and 
avoids the small bias towards a < 1 found in the origi- 
nal modelling. We have demonstrated that bias and red- 
shift space distortions do not change the general shape of 
the correlation function on these scales but can change its 
amplitude. These effects are more apparent in the power 
spectrum approach. The range of wavenumbers used in fit- 
ting models to the power spectrum ratio described in the 
previous paragraph is wider than we use for the correlation 
function, which gives a longer baseline for any gradients to 
become apparent. Also, when forming the correlation func- 
tion from the power spectrum, different scales are mixed 
together, which weakens any gradients in £,{r). The only rel- 
evant impact of redshift space distortions for the case of 
the correlation function is to increase the damping of the 



Sample id 


^ 


52 


fc*/(/iMpc-i) 


Q 




0.0 


1 


0.118 ±0.010 


1.006 ± 0.008 


DM 


0.5 


1 


0.142 ±0.011 


1.002 ± 0.007 




1.0 


1 


0.172 ±0.015 


1.000 ± 0.006 




0.0 


2.796 


0.141 ± 0.043 


0.997 ±0.019 


Halos 1 


0.5 


5.601 


0.141 ±0.044 


1.004 ±0.019 




1.0 


12.04 


0.132 ±0.039 


0.992 ±0.020 




0.0 


3.323 


0.143 ±0.049 


1.009 ± 0.020 


Halos 2 


0.5 


7.306 


0.143 ±0.053 


0.972 ± 0.020 




1.0 


17.20 


0.132 ±0.046 


0.994 ±0.027 




0.0 


6.214 


0.156 ±0.040 


1.002 ± 0.015 


Halos 3 


0.5 


13.29 


0.157 ±0.046 


1.009 ±0.015 




1.0 


31.61 


0.132 ±0.054 


1.008 ± 0.023 



Table 3. The results of applyi ng the general fitting procedure 
described in lAneulo et al.l 1120081 ) to the real-space power spectra 
of the dark matter and the halo samples defined in sec. 14.21 The 
first column gives the label of the sample. The second column 
gives the redshift output. Column 3 gives the mean bias factor 
b^ = (Piih(fc)/^mm(fc)). Columns 4 and 5 give the obtained con- 
straints on fc* and a. All quoted values correspond to the mean 
and 68% c.l. on the respective parameter. 



acoustic peak compared with real space, which results in a 
lower value of the smoothing scale fc* being returned by the 
fitting procedure. The correlation function approach also ig- 
nores the absolute amplitude, but does use the shape of the 
correlation function around the acoustic peak, as well as the 
shape and form of the peak itself. 

We now compare how well the two approaches work 
in practice. As we have already seen in Fig. 1141 a detailed 
comparison of the predictions of Eq. (|14p with the mean cor- 
relation function measured from the ensemble of simulations 
shows small discrepancies. These are the source of the small 
bias towards a < 1 found in § 15.21 On scales larger than the 
acoustic peak the model slightly overestimates the ampli- 
tude of ^(r) which follows more closely the predictions from 
linear perturbation theory. On these scales, Cdw(?") is almost 
entirely determined by the first term of Eq. (|14|) (Crocce & 
Scoccimarro 2008; Smith et al. 2008) . These differences may 
indicate that G(fc), the function that controls the damping of 
the harmonic oscillations, may not be correctly described by 
a Gaussian on these scales. The largest differences between 
the model and the results from the simulations are found 
on scales of r < 90 h~^ Mpc. These discrepancies can be re- 
lated to the change in the shape of the power spectrum due 
to nonlinear evolution. The inclusion of the scale-dependent 
correction factor of Eq. (|16|) can alleviate these differences, 
but it is not enough to correct for the small bias in the con- 
straints on a. This implies that the implementation of this 
modelling when dealing with real observational data may 
lead to a slight bias in the constraints on the dark energy 
equation of state wde. On the other hand the mode l based 
on RPT proposed by ICrocce fc Scoccimarrd ([2008!) yields 
tighter constraints and does not suffer from biases towards 
a < 1. 

Fig. [18] shows the effectiveness of the same model ap- 
plied to P{k). The upper panel shows the mean real-space 
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dark matter power spectrum measured from our ensemble of 
simulations (open points) . The dotted lines indicate the vari- 
ance over the different realizations. Discrepancies between 
the linear theory power spectrum (dot-dashed line) and the 
'dewiggled' power spectrum from Eq. dill) (solid line) are 
evident. Non-linear evolution distorts the shape of P{k) in 
a way that can be correctly described by Eq. (|16|) (dashed 
line) with Q = 13 and A = 1.5. The lower panel of Fig. [18] 
shows the ratio of these power spectra to Pnw(fc). The re- 
sults from the numerical simulations have also been divided 
by the non-linear distortion factor. Although P^^ does a 
good job of reproducing the full shape of the power spec- 
trum and the damping of the oscillations, there are small 
discrepancies that show the limitations of the model. 

In order to assess which of the two approaches is su- 
perior when used as a tool to recover unbiased constraints 
from BAO measurements, we have analysed the mean real- 
space power spectrum from the ensemble of simulations for 
the dark matter and the s ame halo sam p les as defined in 
§ 1321 using the method of lAngulo et al.l l|2008l ). The con- 
strain ts on a and fc ^, obt ained in this way are listed in Ta- 
ble 3. lAngulo et al.l (|2008il defined the mean bias factor as 
6^ — {Ph_h{k) / Pinin{k)) . The values of fc'^ obtained for the dif- 
ferent halo samples analysed are also listed in Table 3, and 
are in quite close agreement with those we quote in Table 1 
for the correlation function. For the case of the dark matter, 
the value of fc* obtained from the power sectrum shows the 
same tendency to increase with redshift as was found in the 
analysis of ^(r). For the halo samples, the trend of k-^ with 
redshift is less evident, perhaps due to the scale dependent 
bias seen in P{k). There is also a tendency towards higher 
values of fc^ with increasing mass which is also suggested in 
the results of § 15.21 The precise preferred values of fc* ob- 
tained in this case are in agreement with the ones obtained 
using the correlation function. 

The values for a recovered from both methods are simi- 
lar, but the allowed region for this parameter is bigger when 
only the information from the oscillations in P{k) is used. 
The extra information from the shape of the correlation 
function helps to improve the constraints on a resulting in a 
smaller allowed region for this parameter. This test implies 
that the simple model for the correlation function described 
in § 15.11 can perform better than the methods in which the 
scale of the acoustic oscillations is extracted from the power 
spectrum. In a more realistic situation, when other cosmo- 
logical parameters are included in the analysis, the extra in- 
formation contained in the shape of the correlation function 
would make this approach seem even more attractive, since 
it can also be used to improve constraints on parameters, 
such as f2m. In the analysis of CMB data, this parameter 
shows a strong degeneracy with wde, and therefore a better 
constraint on the matter density would also result in tighter 
bounds on the dark energy equation of state. 

Finally, the mean values of the stretch parameter re- 
covered from the dark matter power spectrum show a small 
bias towards a > 1 which dec rease with increasing redshifts. 
ICrocce fc Scoccimarrd (|200a ). estimated the mode-coupling 
sh ifts that would affect the constraint on a using the method 
of lAngulo et al.l l|200a ). Their predictions are in good agree- 
ment with our results, suggesting that this bias might be 
attributed to the systematic effects introduced by the mode- 
coupling power spectrum Pmc(fc). This deviation from a = 1 



is smaller than the sampling variance for the volume of the 
individual L-BASICC II realizations. Future surveys like Eu- 
clid or ADEPT will cover volumes that are close to the to- 
tal volume simulated by our ensemble {V — 120 Gpc''/i~'^). 
Rescaling the obtained covariance matrix by a factor 1/50 
(which then corresponds to the covariance of the mean cor- 
relation function from our ensemble) we can gauge the effect 
of this bias for these surveys. In this case, the error associ- 
ated with the value of a at 2 = wiU be 0.008/750 ^ 0.001, 
which then implies a 5 sigma detection of the mode-coupling 
shifts for the dark matter. Instead, the results obtained from 
the correlation function applying the modelling of Eq. p9|l 
are more consistent with a — 1, even after re-scaling the 
error to the full volume of our ensemble. This shows that 
with this modelling, the correlation function is more likely 
to produce unbiased constraints that the power spectrum. 
In principle, the simple ansatz of Eq. (|19[) can be extended 
to the analysis of the power spectrum. Nevertheless, due 
to the fact that in Fourier space, redshift space distortions 
and halo bias do show important scale depe ndence, it will 
be necessary to follow an approach similar to lAngulo et al.l 
(2008), dividing the measured power spectrum by a smooth 
function. This sacrifice of information will always have con- 
sequences for the obtained constraints. 



7 CONCLUSIONS 

In this paper we have addressed two main questions: 1) 
What is the relation between the "acoustic peak" in the 
two-point correlation function and the sound horizon scale? 
and 2) Which is the better statistic to use to constrain the 
dark energy equation of state, the correlation function or the 
power spectrum? Future galaxy surveys will cover volumes 
one or more orders of magnitude larger than existing sur- 
veys. For the first time, systematic effects in the appearance 
of the BAO and in their interpretation will be compara- 
ble to the random or sampling variance errors on the mea- 
surement of the two point statistics. It is therefore essential 
to understand these systematics and to improve theoretical 
modelling of the BAO so that they can be used to provide 
optimal and unbiased estimates of the values of cosmological 
parameters. 

It is a common misconception that the location of the 
acoustic peak in the correlation function is identical to the 
size of the sound h orizon at recombination. Indeed, sev- 
eral recent studies (ICuzik et all l2007l : ISmith et all I2OO8I : 
ICrocce fc Scoccimarrd I2OO81 ) have focused on showing that 
the acoustic peak in the correlation function is shifted and 
distorted relative to the prediction of linear perturbation 
theory. What has not been widely appreciated before is that 
even in linear perturbation theory, the centroid or maxi- 
mum of the peak in ^(r) does not coincide with the sound 
horizon scale at the level of accuracy required to fully ex- 
ploit the constraints that can be derived from forthcoming 
galaxy surveys. The size of the error in the distance scale 
which would result from making this incorrect assumption 
is of the order of 1%, which is already close to the random 
error in the distance scale forecast for ongoing surveys. 

In order to use the correlation function on large scales 
to constrain the values of the cosmological parameters, it 
is therefore necessary to model the form of the correlation 
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function. By using a physically motivated model, a natural 
and consistent connection is made between the cosmological 
parameters and the form of the correlation function. There 
are several effects which have to be taken into considera- 
tion to produce a complete model: the non-linear growth of 
fluctuations, the possible scale dependence of the bias be- 
tween the chosen tracer of the density field and the under- 
lying matter fluctuations and the distortions introduced by 
using distance measurements which are inferred from red- 
shifts, commonly referred to as redshift-space distortions. 
In order to assess the impact of these effects on the correla- 
tion function, we used an ensemble of 50 very large volume 
N-body simulations, the L-BASICC II. These si mulations are 
analog ous to the L-BASICC ensemble used bv lAngulo et al.l 
l|2008h for a similar analysis of the acoustic oscillations in 
the power spectrum. We measured the correlation function 
of the dark matter and of different halo samples in both real 
and redshift-space in each of the realizations in the ensem- 
ble. 

It turns out that bias, general non-linear evolution and 
redshift-space distortions are much simpler to deal with 
in the case of the correlation function than they are for 
the power spectrum. On the scales we consider (60 < 
(r//i~^Mpc) < 180), these effects primarily result in a 
change in the amplitude of the measured correlation func- 
tion and do not alter its shape. By "general" non-linear evo- 
lution, we mean the change in the overall "coarse grained" 
shape of the power spectrum. The crucial effect to include in 
the model is the damping of the higher order harmonic os- 
cillations due to non-linear evolution, as given by Eq. (|lip . 
Whilst we find that a simple damping of the linear per- 
turbation theory power spectrum is able to reproduce the 
correlation function we measure from the simulations, some 
discrepancies remain between the model and the simulation 
results on scales both larger and smaller than the acoustic 
peak. These differences lead to slightly biased constraints on 
the stretch factor a, which in turn will yield small biases in 
the estimates of the dark energy equation of state parameter 
uiDB. The model is only marginally improved on incorpo- 
rating an empirical model for non-linear distortions (which 
could also be used to describe a scale-dependent bias). This 
Q-model, introduced by Cole et al. (2005), is widely used in 
the literature but needs to be refined to extract the full in- 
formation from the correlation function. On the other hand, 
the implementation of the model of Eq. (|19p helps to bring 
the constraints on the stretch parameter in closer agreement 
wi th the value a = 1. This par ametrization was proposed 
bv lCrocce fc Scoccimarrd (|2008l ) and is based on RPT. 

To address the question of which two-point statistic is 
the more powerful for extracting the BAO, we repeated the 
power spectrum analysis carried out by Angulo et al. for sev- 
eral of the samples we used to measure correlation function. 
The analysis of Angulo et al. can be viewed as a conser- 
vative approach (see also iPercival et al.ll2007al la). The mea- 
sured power spectrum is divided by a "wiggle-free" reference 
spectrum which is defined using the measured spectrum it- 
self, without any modelling. No information is used about 
the amplitude of the spectrum or its overall shape. By taking 
a ratio in this way, the BAO are isolated and the impact of 
scale dependent effects in the power spectrum is minimized. 
Similar ratios generated using linear perturbation theory, 
and applying some "de-wiggling" or damping of the BAO 



are compared to the measured power ratio. We compared 
the constraints obtained from the model for ^(r) described 
in Sec. 15.11 a nd the ones obtaine d by applying the fitting 
procedure of lAngulo et al.l (|2008l ) . This simple test shows 
that a method that uses the full large-scale shape of the 
correlation function can provide tighter constraints on wde 
than one in which the scale of the acoustic oscillations is 
extracted from the power spectrum. The extra information 
contained in the shape of the correlation function can also 
be extremely useful to constrain other cosmological param- 
eters, such as r2m. In the analysis of CMB data, f2i„ shows 
a strong degeneracy with woe, and so the use of the corre- 
lation function would result in tighter bounds on the dark 
energy equation of state. 

The constraints on the stretch parameter obtained with 
the power spectrum show a small bias towards a > 1 which 
is consistent with the predictions of Croc ce fc Scoccimarrd 
(2008) based on the effect of the mode-coupling power spec- 
trum. This bias will become much more important for future 
galaxy surveys which will have much smaller random errors. 
Instead, the constraints obtained using Eq. p9|l to model the 
correlation function are much more consistent with a = 1 
even for the large volumes that will be covered by these 
surveys. 

The ability of Eq. (|19|l to improve the constraints ob- 
tained from the correlation function suggests that the recent 
advances in the modelling of the aco ustic oscillations using 
renormalized perturbation theory (ICrocce fc Scoccimarrd 
I2OO3) may provide a basis for a more accurate theoreti- 
cal model of the full shape of the correlation function and 
power spectrum. Nonetheless the problem is aggravated by 
the need to model scale dependent galaxy bias, which may 
have already become the strongest limitation on the use of 
large scale structure info rmation to obtain con straints on 
cosmological parameters (jSanchez fc Cole 1120081 ). An accu- 
rate model capable of describing all sources of distortions 
must be implemented if the determinations of the galaxy cor- 
relation function and power spectrum from future surveys 
are to achieve competitive constraints on the dark energy 
equation of state. It is clear that realistic numerical simu- 
lations which are able to model the growth of structure in 
the dark matter and which can incorporate a physical model 
for galaxy formation will be essential to guide the develop- 
ment of such models. The use of both of these statistics will 
give the complementary information required to improve the 
model of Eq. pip , and will hopefully allow the full poten- 
tial of BAO measurements as probes of the nature of dark 
energy to be reached. 
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